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ABSTRACT 


An  analysis  is  presented  vhlch  can  predict  the  coopressible  laminar  boundary 
layer  characteristics  on  a  blunt  flat  plate  vlth  prescribed  surface  temperature 
and  pressure  gradients.  The  analysis  uses  a  Von  Kazaan  Integral  method  to  solve 
the  compressible  momentum  and,  energy  equations.  Problems  usually  encountered 
with  very  high  pressure  gradients  are  avoided  by  employing  a  new  velocity  profile 
vhich  remains  undistorted.  A  first-order  analysis  of  the  effect  of  entropy  gradi¬ 
ents  in  the  invlscid  flow  on  the  boundary  layer  parameters  is  also  given. 
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I.  INTRODUCTION 


The  use  of  blunt  bodies  to  alleviate  the  severe  heating  problems  associated 
with  atmospheric  re-entry  vehicles  has  generated  a  considerable  amount  of  interest 
in  determining  the  boundary  layer  characteristics  about  blunt  bodies  in  hypersonic 
flow.  During  a  large  portion  of  the  re-entry  trajectory  the  boundary  layer  will 
be  entirely  land and  there  is  some  hope  for  an  analytic  solution,  at  least  in 
the  regime  where  the  Prandtl  boundary  layer  approximations  are  valid.  This  report 
presents  an  attempt  at  solving  the  compressible  laminar  boundary  layer  equations 
for  a  blunt  flat  plate  at  a  zero  angle-of -attack.  The  only  attendant  restrictions 
are  these: 


1.  Prandtl 's  boundary  layer  assumptions  are  valid. 

2.  The  viscosity  is  linearly  dependent  on  temperature. 

3.  The  pressure  distribution  on  the  body  is  available. 

4.  The  wall  temperature  distribution  is  defined. 

3.  The  shock-wave  shape  can  be  predicted. 

Consistent  with  the  above  restrictions,  the  analysis  will  be  capable  of  predicting: 

1.  Boundary  layer  profiles  for  temperature,  enthalpy  and  velocity. 

2.  Local  skin  friction  coefficient. 

3.  Local  heat  transfer  coefficient. 

4.  Total-pressure  gradient  along  the  outer  edge  of  the  boundary  layer. 

5.  Variation  of  the  boundary  layer  thickness  along  the  plate. 

6.  Variation  of  the  momentum,  energy,  and  displacement  thicknesses 
along  the  plate. 

The  literature  contains  many  special  or  restricted  solutions  to  the  laminar  com¬ 
pressible  boundary  layer  equations.  These  solutions  usually  contain  various 
combinations  of  the  following  restrictions: 

1.  Zero  pressure  gradient. 

2.  Zero  surface  temperature  gradient. 

3.  Prandtl  number  equal  to  unity. 

4.  Adiabatic  wall  condition. 

3.  Incompressible  flow. 

6.  Very  specific  pressure  gradients. 


Manuscript  released  by  author  19  October  1962  for  publication  as  an  iSD  Technical 
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7.  linear  viscosity- temperature  relationship. 

8.  Boundary  layer  and  shock  nave  coincide. 

One  of  the  earliest  solutions  of  the  conqjressible  laminar  boundary  layer  equations 
vas  given  by  Chapman  and  Rubesin  (l).*  They  vere  the  first  to  use  the  linear 
viscosity- temperature  relationship 

M  sr  CT 

Ms  T  (l-l) 

to  separate  the  energy  and  momentum  equations.  Unfortunately  their  analysis  does 
not  account  for  the  effects  of  pressure  gradient. 

Basically  the  present  analysis  applies  the  Von  Karman  integral  method  to  ob¬ 
tain  the  solution  of  the  conqsressible  laminar  boundary  layer  equations  for  any 
arbitrary  pressure  and  tenqpepature  distributions.  A  good  reference  to  previous 
works  utilizing  the  integral  approach  for  a  solution  to  the  laminar  boundary  layer 
equations  as  well  as  a  critical  analysis  of  their  short-coming  is  contained  in  the 
introductory  section  of  a  paper  by  Beckwith  (2).  The  analysis  of  Beckwith  has  the 
disadvantage  that  it  cannot  be  applied  to  a  problem  with  a  surface  temperature 
gradient.  The  introductory  section  of  NACA  TN-3157  by  Libby  and  Morduchov  (3) 
gives  another  good  review  of  the  papers  using  the  integral  method  to  solve  the 
compressible  laminar  boundary  layer  equations.  The  analysis  of  NACA  TN-3157  is 
only  valid  for  an  axled,  pressure  gradient  with  a  uniform  wall  temperature  or  a 
zero  pressure  gradient  with  an  arbitrary  welII  temperature  distribution. 

The  analysis  of  this  paper,  in  addition  to  being  valid  for  a  flat  plate  with 
a  wide  range  of  pressure  and  temperature  distributions,  can  be  extended  to  encompass 
the  solution  of  three-dimensional  axially- symmetric  blunt  bodies  at  zero  angle-of- 
attack,  as  well  as  blunt  (or  sharp)  flat  plates  at  moderate  angles-of -attack. 

The  only  severe  restriction  on  the  solution  is  the  linear  dependence  of  the 
viscosity  on  temperature.  The  analysis  by  Chapman  and  Rubesin  (l)  to  define  the 
error  involved  by  using  the  linear  relationship  instead  of  the  more  accurate 
Sutherland  law  concluded  that  "regardless  of  Mach  number,  free-stream  temperature, 
or  average  surface  temperature,  the  approximation 

^5  Cl 

Moo  To 

can  be  made  to  yield  nearly  correct  values  for  the  viscosity  in  the  important 
region  near  the  surface."  They  then  show  that  the  approximate  viscosity  law  does 
not  introduce  an  error  of  more  than  5  or  6  per  cent  over  the  range  of  free  stream 
temperatures  from  72°R  to  6i+8°R  and  a  Mach  number  over  the  plate  of  5<  The  maximum 
error  in  the  boundary  layer  parameters  obtained  by  the  analysis  of  this  paper  attri¬ 
buted  to  using  the  linear  viscosity-temperature  relationship  insteeul  of  the  Suther¬ 
land  equation  could  be  expected  then  to  be  around  5  or  6  per  cent  since  in  hyper¬ 
sonic  flow  the  Mach  number  over  the  plate  never  gets  much  higher  than  5  and  there 
is  reason  to  believe  that  the  error  is  not  a  strong  function  of  temperature. 


^Numbers  in  parentheses  denote  references  listed  in  the  bibliography. 
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To  start  the  analysis  the  shock- vave  shape,  the  surface  temperature  distri¬ 
bution,  and  the  pressure  distribution  existing  on  the  body  must  be  known.  Part 
of  the  problem  of  solving  for  the  boundary  layer  characteristics  is  to  define 
both  the  pressure  distribution  and  shock-vave  shape  consistent  with  the  body 
shape  and  free- stream  Mach  number.  One  method  of  obtaining  the  pressure  distri¬ 
bution  about  a  blunt  flat  plate  at  very  high  free- stream  Mach  numbers  is  taken  up 
in  Section  II.  A  method  of  predicting  the  shoe*- vave  shape  is  discussed  in 
Section  III.  The  methods  selected  for  defining  the  shock-vave  shape  and  pressure 
distribution  vill  be  pointed  towards  applications  in  hypersonic  flow  vhere  at 
present  the  greatest  interest  seems  to  lie.  This  does  not  in  any  vay  detract  from 
the  usefulness  of  the  analysis  at  any  steady-state  supersonic  flow  condition. 


II.  ANALYTIC  METHOD  FOR  PREDICTING  THE  PRESSURE  DISTRIBUTION 
ABOUT  BLUNT  BODIES  AT  ZERO  ANGLE- OF- ATTACK 

The  pressure  distribution  over  the  body  must  be  known  before  the  boundary 
layer  equations  can  be  solved  by  the  Von  Karman  integral  method.  In  the  past  an 
integral  analysis  itself,  to  which  another  constraint  has  been  added,  has  been 
used  to  predict  the  pressure  distribution  (4).  This  added  constraint  vas  to 
assume  that  the  shock  vave  lies  right  on  the  boundary  layer.  This  constraint  has 
been  criticized  by  some  (5)  as  not  being  representative  of  the  physical  conditions, 
consequently  it  is  not  used  in  this  analysis.  Instead,  a  method  of  predicting  the 
pressure  distribution  will  be  used  which  has  been  found  to  be  very  successful  (6). 
The  method  basically  consists  of  predicting  the  pressure  distribution  on  the  cylin¬ 
drical  leading  edge  using  the  Prandtl-Meyer  plus  Newtonian  analysis  as  proposed  by 
Lees  and  Kubota  (6)  and  using  the  hypersonic  blast-vave  analogy  to  define  the 
pressure  distribution  on  the  afterbody. 

In  the  nose  region  it  has  been  found  experimentally  that  the  Newtonian  corpus¬ 
cular  theory  as  modified  by  Lees  adequately  predicts  the  pressure  distribution  on 
a  cylindrical  leading  edge  up  to  a  point  about  55°  from  the  shoulder.  Lees  and 
Kubota  (6)  in  attempting  to  predict  the  pressure  distribution  to  the  shoulder  tried 
matching  a  Prandtl-Meyer  expansion  to  the  Newtonian  solution.  The  point  at  which 
the  two  methods  were  Joined  vas  defined  to  be  that  point  at  which  the  pressures  and 
predicted  by  both  analyses  are  the  same.  The  above  conditions  together  with 
the  condition  that  the  total  pressure  was  assumed  constant  around  to  the  shoulder 
allowed  one  to  calculate  the  angle  and  Mach  number  Mm?  at  the  match  point  as 

a  function  of  the  free- stream  Mach  number.  A  plot  of  the  variation  in  and  Map 
with  free- stream  Mach  number  is  found  in  reference  (7) .  Fortunately  the  pressure 
distribution  generated  by  such  an  analysis  has  been  shown  experimentally  to  predict 
the  inviscid  pressures  about  a  blunt  body  around  to  the  shoulder.  It  will  not 
account  for  any  induced  pressure  effects. 

A  comparison  of  the  non-dimensional  pressure  distribution,  f*/P,  as  predicted 
by  the  Prandtl-Meyer  plus  modified  Newtonian  analysis,  is  made  in  Figures  1  and  2. 
The  experimental  data  at  M^  -  6.86  (Figure  l)  were  obtained  from  reference  (8)  and 
the  experimental  data  at  -  12.28  (Figure  2)  from  reference  (9).  From  Figures  1 
and  2  it  can  be  seen  that  the  modified  Newtonian  plus  Prandtl-Meyer  analysis  can  be 
considered  as  a  satisfactory  approximation  of  the  pressure  distribution  around  the 
cylindrical  leading  edge.  It  might  be  noted  that  in  both  cases,  the  blunt  body  vas 
a  complete  cylinder  and  the  effect  of  an  afterbody  is  not  well  known.  In  reference 
(9),  a  cylindrically  blunted  flat  plate  vas  tested  at  moderate  angxes-of- attack; 
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Figure  3  presents  s  portion  of  the  data  showing  that  the  forebody  pressures  are  not 
Influenced  by  these  angles- of- attack.  Somewhat  conflicting  results  were  obtained 
In  reference  (10),  a  series  of  tests  on  a  flat  plate  of  slightly  different  config¬ 
uration  than  the  blunted  slab  of  reference  (9).  As  shown  in  Figure  3>  it  appears 
that  the  distributions  are  affected  by  the  afterbody  at  angles- of- attack.  Whether 
these  differences  can  be  attributed  to  experimental  scatter,  Reynolds  number,  or 
to  model  configuration,  or  whether  the  forebody  pressures  are  indeed  influenced  by 
afterbody  shape  will  have  to  be  resolved  by  further  experiment. 

It  must  be  kept  in  mind,  however,  that  this  representation  neglects  the  effects 
of  Reynolds  number  as  well  as  surface  temperature  distribution.  Extreme  care  must 
be  taken  to  insure  that  a  method  of  predicting  pressures  around  blunt  leading  edges 
is  not  used  in  situations  where  it  is  not  applicable,  i.e.,  flow  conditions  where 
severe  merging  of  the  boundary  layer  and  shock  wave  caused  an  induced  pressure  to 
alter  the  inviscid  (Prandtl-Meyer  plus  Newtonian)  solution  appreciably.  Along  this 
line  the  analysis  of  reference  (ll)  will  prove  invaluable. 

The  pressures  on  the  afterbody  can  be  determined  from  the  hypersonic  flow-blast 
wave  analogy.  It  has  been  demonstrated  (12)  that  the  equations  of  motion  of  the 
hypersonic  «m*n  disturbance  theory  (13)  are  analogous  to  the  equations  of  motion 
for  unsteady  shock  wave  flow  (l4,  15,  16,  17,  18) .  The  analysis  of  Lees  and  Kubota 
(6)  indicates  that  in  particular  the  analogy  is  a  valid  one  for  the  case  of  hyper¬ 
sonic  flow  over  a  flat  plate  with  a  cylindrical  leading  edge.  The  analogy  is  an  im¬ 
portant  one  when  it  is  recognized  that  the  blast- wave  equations  (equations  of  motion 
for  unsteady  shock  waves)  have  been  solved  many  times  (l4,  15,  l6,  17,  18)  and  as  a 
consequence  of  the  analogy  these  solutions  represent  a  solution  to  the  equation  of 
motion  of  the  hypersonic  small  disturbance  theory.  Lukasiewicz  (19)  develops  equa¬ 
tions  for  defining  the  surface  pressures  and  shock-wave  shape  applicable  to  a  blunt 
flat  plate  from  the  solution  of  the  blast-wave  equations  of  Sakurai  (14,  15).  Die 
paper  by  Lukasiewicz  is  also  an  excellent  survey  paper  pertaining  to  the  hypersonic 
flow-blast  wave  analogy. 

Sakurai  presents  two  solutions  to  the  blast- wave  equations.  In  both  solutions 
the  basic  property  of  the  analysis  is  the  power  series  representation  for  pt  and 
U  in  terms  of  the  inverse  shock  wave  Mach  number  £*-//*  wj  •  For  very  strong  ex¬ 
plosions  the  speed  of  sound,  CL,  will  be  very  much  smaller  than  the  wave  speed,  U.^ 

Die  first  solution  keeps  the  linear  terms  in  {‘lJv.J]  (l4)  and  the  second  solution 
the  terms  up  to  „-Ja  (15)  •  In  addition  to  the  approximate  solutions  of 
Sakurai  (l4,  15)  the  blast-wave  equations  have  an  exact  solution  which  was  deter¬ 
mined  by  Sedov  (l6) .  Die  exact  solution  of  Sedov  agrees  with  Sakurai ' s  first-order 
solution,  being  5.2  per  cent  larger  at  all  values  of  x/d,  but  does  not  agree  nearly 
as  well  with  Sakurai 's  second-order  solution.  A  comparison  of  the  solutions  of 
Sedov  (16)  and  Salurai  (l4,  15)  with  experimental  data  is  made  in  Figure  4.  It  is 
evident  from  Figure  4  that  the  exact  solution  of  Sedov  does  not  represent  the  data 
as  well  as  Sakurai 's  second-order  solution.  Diis  is  not  too  alarming,  however,  since 
in  developing  the  hypersonic  small  disturbance  theory  the  linearization  process 
could  easily  account  for  the  discrepancy  between  the  experimental  data  and  the  exact 
solution  of  the  linearized  equations  of  motion.  It  might  also  be  noted  that  the 
over-dl  agreement  i6  the  poorest  in  the  vicinity  of  the  nose  region.  This  is  what 
might  be  expected  because  it  is  in  the  nose  region  that  the  assumption  that  the 
body  slopes  be  small,  which  effectively  means  that  the  velocity  disturbances  are 
small,  is  invalid.  All  available  correlation  of  theory  with  experimental  data  in¬ 
dicates  that  the  blast-wave  theory  is  valid  approximately  for  X*  Z-  2  or  3. 
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The  equation  which  best  defines  the  pressure  distribution  over  a  blunt  flat 
plate,  as  obtained  from  Sakural's  second-order  solution,  is 

^=12,M-t%]S+OS6  (2-1) 

The  term  CoM in  equation  (2-l)  is  analogous  to  the  initial  energy  of  the  explosion 
in  the  blast-wave  equations  arai  its  full  significance  is  explained  in  reference  (6). 
The  above  theory  as  applied  to  hypersonic  flow  has  the  disadvantage  that  it  does 
not  account  for  any  effects  of  heat  transfer  or  surface  temperature  distribution. 
Fortunately  it  seems  from  an  analysis  of  experimental  data  that  this  is  not  too 
severe  a  limitation,  at  least  not  of  the  order  of  magnitude  of  the  next  problem. 

The  initial  hypersonic  analysis  assumed  inviscid  flow  and  at  low  free-stream  Rey¬ 
nolds  numbers  it  might  be  expected  that  an  interaction  between  the  boundary  layer 
and  shock  wave  could  alter  the  inviscid  pressure  distribution  over  the  plate.  This 
Reynolds  number  effect  is  illustrated  in  Figure  5 >  where  experimental  pressure  data 
from  reference  (9)  have  been  plotted.  It  can  be  seen  that  the  experimental  data  at 
a  Reynolds  number  of  11,900  do  not  agree  with  the  second-order  blast-wave  theory  as 
well  as  the  data  taken  at  the  higher  Reynolds  numbers. 

The  pressure  distribution  in  the  vicinity  of  the  shoulder,  where  both  methods 
discussed  are  not  applicable,  is  obtained  by  joining  the  data  from  the  Prandtl-Meyer 
expansion  with  that  of  the  blast-wave  theory  with  a  reasonably  smooth  curve.  It  is 
noteworthy  to  mention  that  if  the  data  for  the  pressure  over  the  plate  from  both 
theoretical  analyses  are  expressed  in  the  form  v*  X*  ,  they  will  generate  one 

curve  for  NA,e8.  This  would  lead  one  to  the  conclfision  that  the  predominant  para¬ 
meter  defining  the  pressure  on  the  blunt  flat  plate  is  the  impact  pressure  down¬ 
stream  of  the  normal  shock  and  not  the  free-stream  Mach  number.  This  conclusion 
is  borne  out  for  the  leading  edge  region  by  the  experimental  data  as  shown  in  Figure 
6.  The  pressures  on  the  afterbody  do  not  correlate  in  quite  the  same  manner  as 
shown  in  Figure  7  where  the  data  of  reference  (9)  are  plotted  for  three  different 
Mach  numbers  all  with  the  same  and  ~T^.  .  The  correlation  is  not  as  Mach-number 

independent  as  the  theory  would  indicate.* 

In  addition  to  the  pressure  distribution  on  the  body  the  pressure  gradient 
term  must  also  be  known  in  order  to  obtain  a  solution.  In  analyzing  the 

gradient  term  from  the  Prandtl-Meyer  plus  Newtonian  analysis  it  was  noticed  that, 
although  the  pressure  gradients  at  the  match  point  are  equal,  the  curvature  or 
second  derivative  of  the  pressure  with  respect  to  x  is  not  continuous.  In  the  non- 
dimensional  plot  of  the  pressure  gradient  versus  X*  (Figure  8)  this  appears 

as  a  discontinuity  in  the  slope.  Figures  8  and  9  present  a  comparison  between  ex¬ 
perimentally  and  theoretically  determined  pressure  gradients.  From  a  comparison  of 
the  data  presented  in  Figures  8  and  9  it  can  be  concluded  that  for  the  purposes  of 
a  boundary  layer  analysis  the  theoretically  determined  pressure  gradients  are  an 
acceptable  approximation  to  the  experimental  data. 


III.  SELECTION  OF  THE  METHOD  FOR  PREDICTING 
THE  SHOCK-WAVE  SHAPE 

As  previously  mentioned,  in  order  to  arrive  at  a  numerical  solution  for  the 
boundary  layer  characteristics  the  shape  of  the  shock  wave  in  the  vicinity  qf  the 
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stagnation  point  mist  be  known,  The  stagnation  region  Is  of  primary  Interest 
because  all  the  mass  flow  in  the  boundary  layer  will  have  passed  through  the 
strong  portion  of  the  shock  wave,  that  portion  of  the  shock  wave  which  is  in  the 
vicinity  of  the  forward  stagnation  point.  Three  methods  were  examined  in  detail 
before  one  was  selected  to  predict  the  shock-wave  shape  around  the  test  body 
(20,  21,  22).  A  plot  of  the  shock  wave  in  the  vicinity  of  the  stagnation  point 
as  predicted  by  the  methods  of  references  20,  21,  and  22  is  shown  in  Figure  10. 
Fortunately  the  three  methods,  although  differing  in  complexity,  predict  basically 
the  same  shock-wave  shape  near  the  stagnation  point.  Consequently,  since  there 
was  no  appreciable  difference  among  the  three  methods,  the  one  easiest  to  adapt 
to  our  particular  situation  was  selected.  Tills  was  the  method  of  reference  (20). 
The  shock-wave  shape  is  represented  by 

Z*=  2X*  ['  +3oTrtl  ,  (3-1) 

where  X  is  the  normal  shock  density  ratio. 

References  20,  21,  and  22  are  by  no  means  the  only  papers  available  for 
predicting  shock-wave  shapes.  These  references  were  selected  because  each  satis¬ 
fied  the  high  Mach  number  restriction.  If  an  analysis  is  to  be  carried  out  at  a 
lover  Mach  number,  the  particular  method  chosen  to  predict  the  shape  at  Mach  12 
cannot  be  applied  since  it  is  an  hypersonic  approximation.  But  there  are  other 
methods  for  predicting  shock-wave  shape  (23)  which  vould  work  equally  well  with 
this  analysis. 


IV.  THEORETICAL  ANALYSIS 

This  analysis  presents  a  method  of  solving  the  compressible  boundary  layer 
equations  for  a  blunt  flat  plate  at  zero  angle-of-attack.  Hie  main  purpose  of 
the  analysis  is  to  determine  the  boundary  layer  characteristics  on  a  blunt  flat 
plate,  without  referring  to  any  experimental  data  if  possible.  Hie  only  restric¬ 
tions  imposed  on  the  analysis  are  the  following: 

1.  The  pressure  distribution  on  the  body  1b  available. 

2.  The  shock-wave  shape  can  be  predicted. 

3.  The  vail  temperature  distribution  is  defined. 

U.  The  viscosity  is  linearly  dependent  on  temperature. 

5.  Prandti's  boundary  layer  assumptions  are  valid. 

The  restrictions  (l)  suid  (2)  were  anticipated  in  sections  II  and  III  where  ana¬ 
lytical  methods  for  obtaining  the  shock-wave  shape  and  the  pressure  distribution 
on  the  body  were  discussed.  Restriction  (3)  can  toe  relaxed  if  the  adiabatic  wall 
condition  is  analyzed  and  in  effect  the  adiabatic  wall  temperatures  can  be  defined 
for  a  flat  plate  with  a  pressure  gradient.  Restriction  (4)  was  discussed  in  Section 
I  and  it  was  shown  that  the  error  to  be  expected  is  less  than  6  per  cent.  Inherent 
in  (5)  is  the  fact  that  =  0.  This  expression  for  fails  in  the  region 
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vbere  the  boundary  layer  thickness  Is  large  compared,  to  the  local  radius  of  curvature 
of  the  body,  lhe  cylindrical  leading  edge  of  the  blunt  flat  plate  under  analysis 
is  highly  curved,  and  consequently  the  analysis  to  follow  will  be  valid  only  for 
conditions  vbere  -ct  |  . 

Consider  the  flow  model  shown  in  Figure  11.  All  the  mass  in  the  boundary  layer 
at  the  station  Xrt  has  passed  through  the  shock  wave  between  the  stagnation  stream¬ 
line  and  the  streamline  .  Since  the  streamline  is  being  followed  from  the 
point  where  it  crosses  the  shock  wave  to  the  point  where  it  enters  the  boundary 
layer,  the  mass  in  the  streamtube  crossing  the  shock  wave  between  and  must 
be  equal  to  the  mass  in  the  boundary  layer  at  Kr\ .  If  the  mass  crossing  the  shock 
wave  between  and  is 


m.-r  p.u,£t 

and  the  mass  in  the  boundary  layer  at  is 


(4-1) 


(4-2) 


then  by  equating  the 
obtained. 


vhere  the  subscript  refers  to  conditions  at  the  outer  edge  of  the  boundary 
layer.  The  term  Ps^s  can  be  taken  under  the  integral  sign  because  it  is  not  a 
'y' -dependent  term. 


two  expressions  for  mass  flow  the  following  equation  is 

.  ?.**  [*  -£“  Am  . 

*"  j>,U.  J0  (4-5) 


In  order  to  make  the  analysis  effectively  one  of  constant  density,  the 
Dorodnitzyn  transformation  is  used.  It  can  be  written. 


Pj>s  =  An.  (4-4) 

or  if  v\_»  ?  ,  where  A  is  the  boundary  layer  thickness  in  the  transformed  co¬ 

ordinate  sysras,  it  can  be  written  as 


(4-5) 


Transforming  the  co-ordinates  in  equation  (4-5)  in  accordance  with  equation  (4-5) 
the  following  equation  results: 


(4-6) 


Nov  following  the  method  of  Von  Kerman  the  velocity  ratio  —  is  assumed  to  be 
expressible  as  Mi 
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(4-7) 


—  *  4  a 

Mg 

In  effect  the  unknown  boundary  layer  profile  baa  been  traded  for  five  unknown  con¬ 
stants  in  a  fourth-order  polynominal  in  .  To  solve  for  the  five  coefficients 
five  boundary  conditions  on  the  profile  are  needed.  It  has  been  shown  (3)  that 
it  is  advantageous  to  distribute  the  boundary  conditions  equally  between  the  wall 
and  the  outer  edge.  For  this  analysis  the  boundary  conditions  will  be  chosen  as 


1. 

t^-<5  , 

u  =  x 

2. 

Cr\m  0  , 

6%i 

3. 

^-3  Cxx-O  , 

4. 

U  =  o 

(4-8) 

The  fifth  boundary  condition  is  obtained  by  evaluating  the  momentum  equation  at 
the  wall.  If  this  is  done,  one  obtains  the  following  expression: 


The  first  three  boundary  conditions  of  equation  (4-8)  state  that  the  velocity  at 
the  outer  edge  of  the  boundary  layer  is  Mj  and  that  there  is  no  change  in  the 
first  or  second  derivative  of  u  with  respect  to  y  at  the  outer  edge  of  the  bound¬ 
ary  layer.  When  there  are  entropy  gradients  in  the  inviscid  flow  field,  the  second 
and  third  boundary  condition  will  have  to  be  modified.  This  will  be  further  dis¬ 
cussed  in  Section  VII.  Hie  fourth  boundary  condition  states  that  the  velocity 
adjacent  to  the  wall  is  identically  zero. 

As  is  usually  done,  the  boundary  conditions  of  equation  (4-8)  are  used  to 
define  coefficients  I  ,°^3  and  as  a  function  of  and  then  equation 

(4-9)  i3  considered  the  defining  equation  for  .  If  this  is  done,  the  follow¬ 
ing  polynominal  for  %x  is  obtained, 

© 


an:ar£*rC  ♦  ■^[irvav3a-,i}  t  ^ 

where  -  2  Introducing  the  expression  for  of  equation  (4-10)  in 

equation  (4-o)  and  evaluating  the  resulting  integral,  the  following  equation  is 
obtained, 


A 

d 


(4-11) 
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In  evaluating  the  integral  it  has  been  assumed  that  .A*  +<.*0 . 

This  will  now  he  shown  to  be  the  case  by  evaluating  the  fifth  boundary  condition, 
equation  (4-9).  It  is  at  this  point  in  the  analysis  that  the  Chapman- Rube sin 
viscosity-temperature  relationship  is  used. 

The  most  accurate  viscosity  relationship  with  which  the  term  ( x4T.X-*could 
be  determined  is  that  of  Sutherland,  at  least  in  the  range  of  temperatvfl*e6  likely 
to  be  encountered  on  a  blunt  body  in  hypersonic  flow;  however,  it  is  difficult  to 
incorporate  into  the  integral  analysis.  The  problems  presented  by  the  Sutherland 
law  can  be  circumvented  by  using  the  Chapman-Rubesin  viscosity- temperature  relation¬ 
ship.  It  is  expected  that  their  relationship  will  introduce  errors  no  larger  than 
5  or  6  per  cent. 


Using  the  result  that  “fyoy  *0  (obtained  from  an  order  of  magnitude  analysis 
of  the  y-momentum  equation),  theflDorodnitzyn  transformation,  equation  (4-5),  can  be 
written 


1, 


A 


Utilizing  the  Chapman-Rubesin  viscosity  model. 


(4-12) 


Ms 


=•  C 


(4-13) 


in  conjunction  with  the  Dorodnitzyn  transformation,  equation  (4-12),  and  the  poly- 
nominal  expression  for  the  velocity,  equation  (4-10),  the  momentum  equation  at 
the  wall,  equation  (4-9),  can  be  used  to  definesA  in  the  following  manner, 


■A.* 


X.  s 

Us/*-* 


(4-14) 


It  is  evident  from  the  above  equation  that  A  ,  the  velocity  boundary  layer  shape 
parameter,  is  only  a  function  of  x  and  not  y  and  consequently  the  integration  per¬ 
formed  to  obtain  equation  (4-11)  is  valid.  In  many  applications  equation  (4-l4) 
will  be  used  as  an  expression  for  A*  as  a  function  of  x  in  the  following  form. 


A  U, T„* 

C-  <*%*) 


(4-15) 


The  compressible  momentum  integral  equation  can  be  written,  following  the 
analysis  of  Schlichting  (24),  as 

^--0,3x1  ©  'J  (w« 

where  d  is  the  momentum  thickness  and  £  the  displacement  thickness.  Evaluating 
fr0B  e(^ua't'ion  (4-10)  with  the  help  of  the  Dorodnitzyn  transformation, 
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equation  (4-12),  using  the  compressible  Bernoulli  equation, 

4!±=_1.  /.dpi 

3x  fit^s  L  dx  J  ,  (it-iy) 


to  eliminate  the  velocity  gradient  and  then  multiplying  both  sides  of  equation 
(4-16)  by  the  momentum  thickness,  Q  ,  results  in  the  following  equation: 


i£l=  J  !(£*)!  - 

6X  A  L  T*  a  J 


(4-18) 


Equation  (4-18)  is  the  compressible  momentum  integral  equation  in  the  form  it  will 
be  used  in  the  numerical  calculations. 

The  terms  and  which  appear  in  equation  (4-18)  have  yet  to  be  evalu¬ 

ated.  Then  momentum  thickness,  O  ,  can  be  defined  as 


(4-19) 


Applying  the  Dorodnitzyn  transformation  in  the  form  of  equation  (4-5)  and  then 
evaluating  the  resulting  integral,  one  obtains  the  following  expression  for  the 
momentum  thickness: 


& 

A 


3«S 


L 


37 


.  A. 


*] 

i  J  > 


(4-20) 


where  it  is  noted  that  -j  .has  the  same  form  as  in  incompressible  flow.  If  the 
term  <5*/©  is  written  as  /°/A  ,  then  it  is  evident  that  an  expression  for 
the  non-dimensional  displacement  thickness  sX  is  still  needed. 


% 


The  evaluation  of  %  is  not  as  straightforward  as  was  the  evaluation  of 
.  The  displacement  thickness  in  the  rjj-plane  can  be  written  as 


i 

A 


m-sjH 


where  again  the  fact  that  = 

already  been  shown  to  equal  w  * 

vto  1 


may  be  written, 


L'o 


(4-21) 

has  been  used.  The  integral  has 

and  with  this  the  displacement  thickness 


k  -n 

c  Jo  ' 


7  A. 

<-  *  IO~  V£© 


(4-22) 


This  is  the  expression  that  will  be  used  to  evaluate  for  use  in  the  numerical 

calculations.  More  will  be  said  later  concerning  the  evaluation  of  the  integral 

S'.?  • 
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The  term  &  which  appears  in  the  momentum  integral  equation  is  evaluated 
from  the  expressions  for  (equation  (4-20))  and  A*  (equation  (4-15))  in  the 
following  manner, 


(4-25) 


Again  consider  the  sketch  shown  in  Figure  11.  Since  the  total  enthalpy  of 
the  flow  is  unaffected  by  passing  through  the  shock  wave,  the  streamline 
enters  the  boundary  layer  at  with  free- stream  total  enthalpy.  In  contrast 
to  the  conservation  of  mass  analysis,  where  the  free-stream  mass  and  the  mass 
at  Xn  were  equated,  the  energy  in  the  free  stream  between  V,  and  cannot  be 
equated  to  the  energy  in  the  boundary  layer  at  because  there  has  been  a  net 
heat  transfer  to  the  wall.  It  can  be  accounted  for  by  summing  up  the  energy  lost 
to  the  wall  between  X*0  and  The  energy  balance  can  then  be  written 


(4-24) 


where  p.  U.V\<k2,  is  the  free-stream  energy  in  the  unit  streamtube  bounded  by 
y  and  y  .  C* is  the  energy  contained  in  the  boundary  layer  at  station 
and  JC*  *****  ^  the  heat  lost  to  the  wall  (or  added  to  the  boundary  layer 
depending  on  the "wall  temperatures).  If  both  sides  of  equation  (4-24)  are  divided 
by  and  then  the  Dorodnitzyn  transformation  (equation  (4-4))  is  applied, 

the  resulting  equation  is 


(4-25) 


where  equation  (4-6)  has  been  used  to  simplify  the  above  expression. 


If  the  energy  thickness  is  defined  as 


e-*  f/EU'-fe]  h- 


(4-26) 


then  can  be  obtained,  by  applying  the  transformation  from  4  ton  (equation 

(4-4)),  as  w 

!“'  i  “r  -K.  1  ^  • 


(4-27) 


This  definition  of  £  allows  the  energy  equation  (equation  (4-25))  to  be  written 
as 


(4-28) 
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The  parameter  Z  will  be  evaluated  after  a  suitable  polynominal  expression 
for  '''vVh,  has  been  defined.  As  it  stands,  equation  (4-28)  is  basically  the 
energy  integral  equation  and  will  be  used  to  define  the  heat  transfer  character¬ 
istics  of  the  body. 

For  purposes  of  evaluating  the  parameter  »•  the  total  enthalpy  in  the 
boundary  layer  can  be  assumed  to  have  the  form 


where 

Equation  (4-30)  represents  the  total  enthalpy  distribution  for  the  situation  where 
the  wall  is  at  the  adiabatic  temperature.  The  expression  for  the  total  enthalpy 
profile  through  a  boundary  layer  with  heat  transfer  is  represented  as  perturbation 
of  the  adiabatic  wall  condition  in  equation  (4-29). 

The  function  fecri)  is  assumed  to  have  the  form 


(4-29) 


(4-30) 


fee 


~ 


i«o 


The  following  boundary  conditions  define  the  coefficients  in  equation  (4-31): 


1. 

Ctvm) 

Vi*  » w 
©  1 

2. 

^•8  Or\,»0 

’  j"*  —  o 
cVtv. 

5. 

4-°  . 

r  ^«ur  . 

4. 

V 

4-0  C<\»o) 

0 

ii 

5  fa 

(4-32) 

Physically  the  first  boundary  condition  expresses  the  fact  that  at  the  outer 
edge  of  the  boundary  layer  the  total  enthalpy  of  the  equilibrium  boundary  layer  is 
equal  to  the  free-stream  total  enthalpy.  The  secondary  boundary  condition  is  tied 
into  the  first  in  that  it  says  all  of  the  fluid  beyond  the  outer  edge  of  the  boundary 
layer  has  a  constant  total  enthalpy  equal  to  the  free-stream  total  enthalpy.  Hie 
third  boundary  condition  states  that  equilibrium  total  enthalpy  of  the  fluid  ad¬ 
jacent  to  the  wall  is  obtained  by  considering  the  fluid  at  the  wall  temperature  and 
pressure.  The  fourth  boundary  condition  arises  from  the  fact  that  for  zero  heat 
transfer  to  the  wall  the  temperature  gradient  at  the  wall  ar/.  «c>  ,  and  con¬ 
sequently  /L_  m  O  • 
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Applying  the  boundary  conditions  of  equation  (4-52)  to  the  expression  for  the 
equilibrium  total  enthalpy  distribution  the  following  equation  is  obtained  for  V\^ 


a  ♦  [.Vw  -  ♦  Si  ,  (4-55) 

The  term  arises  because  there  were  only  four  boundary  conditions  specified  to 
determine  five  coefficients,  and  just  as  was  done  with  the  velocity  profile  four 
of  the  coefficients  were  determined  as  a  function  of  the  fifth. 

Again  a  fourth-order  polynominal  of  the  form 


c^c/O 


&  ^ 


(4-3-0 


is  used  to  evaluate  the  distribution  of  total  enthalpy  through  the  boundary  layer. 
The  boundary  conditions  which  define  the  coefficients  in  equation  (4-54)  are: 


1. 

,  S-Vu,  • 

2. 

rs 

C  t  ”3 

_  o  . 

3. 

Vn-Vw  .  (4-55) 

The  first  boundary  condition  stems  from  the  fact  that  at  the  outer  edge  of  the 
boundary  layer  the  total  enthalpy  of  the  boundary  layer  reaches  the  free -stream 
value.  The  second  condition  indicates  that  the  local  rate  of  change  of  the  total 
enthalpy  is  zero  at  the  outer  edge  of  the  boundary  layer.  The  third  condition 
states  that  at  the  wall  the  total  enthalpy  of  the  fluid  may  be  obtained  from  the 
surface  pressure  and  wall  temperature.  This  time  there  are  only  three  boundary 
conditions  with  which  to  solve  for  five  unknown  coefficients.  Consequently  any 
three  of  the  coefficients  can  be  determined  as  a  function  of  the  other  two.  Carry¬ 
ing  out  the  indicated  calculation  the  expression  for  V\^C'j5  becomes, 


+**<-*t>  ♦  5  CVwOl  , 

*  (4-56) 

where  equation  (4-35)  has  been  substituted  for  .  There  are  three  coef¬ 
ficients  in  equation  (4-36)  which  still  must  be  determined, ^  ,  £  ,  and  . 

The  expressions  for  and  Sx  can  obtained  by  considering  the  energy 
equation  evaluated  at  the  mil, 


1? 


(4-37) 


Expanding  equation  (4-37),  using  the  Chapman-Rube sin  expression  to  represent  the 
thermal  conductivity,  k,  and  applying  the  Dorodnitzyn  transformation  (equation 
(4-12)),  the  expression  that  will  be  used  to  determine  and  £  becomes 


\tl]  + 

t 

—  1  =  o 

1 

L**]„ 

(4-38) 

At  this  point  in  the  analysis  the  above  equation  will  be  used  to  determine  ^  . 
Using  the  basic  definition  for  the  specific  heat  at  constant  pressure, Cy  — 
to  convert  from  a  temperature  derivative  to  a  enthalpy  derivative  and  then  using 
the  equilibrium  enthalpy  distribution  (equation  (4-33))  to  evaluate  the  enthalpy 
derivative  at  the  vail  the  expression  for  $  becomes 

where  was  evaluated .  using  the  velocity  distribution  function  (equation 

(4-10)).  Applying  the  same  procedure  but  using  the  enthalpy  distribution  of  equa¬ 
tion  (4-36)  to  evaluate  the  enthalpy  derivative  the  expression  for  becomes 


r.  -  (§L  t 


*  V^flu*T  -  Vu 

zcf 


(4-40) 


f. 


(4-41) 


then  the  total  enthalpy  distribution  through  the  boundary  layer  may  be  expressed 
as  a  function  of  the  coefficient  to*  and  in  the  following  manner, 

«  <c ♦  K,Ca'  V>„  V><.( 


W  {(.'-% } 


(4-42) 


The  functions  and  are  the  only  terms  dependent  on  the 

variable  v^_ in  the  expression  for  the  total  enthalpy  distribution.  If  equation 
(4-42)  is  divided  by  Vv*  and  the  resulting  expression  substituted  for  in 
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©u 

equation  ( 4-27 ) >  then  the  expression  for  can  be  developed  to  the  point  where 

it  can  be  written  as 

ft,,  ,iAi 

A  V\^  J4£0  1*4  1  “i  I  Ztto  *>oZ4  J  ~  V.  &T Su 

-  bsAf  . 

Vv,  l 

The  expression  for  the  shock-wave  shape,  which  was  developed  in  Section  III, 
will  now  be  used  to  obtain  some  relationships  which  are  necessary  to  complete  a 
numerical  solution.  The  equation  defining  the  shock  wave,  equation  (3-1),  can  be 
differentiated  once  to  obtain  the  slope. 


(4-43) 


cU*  _  TAVi  Q  _  1 

ax,  M  ,  (4-mo 

where  @ur  is  defined  in  Figure  11.  Equation  (4-44)  can  be  considered  as  a  relation¬ 
ship  between  dia  and  2,  ,  wit);  another  relationship  between  61*  and  existing 
in  the  oblique  shock  tables  (25),  or  effectively  there  exists  a  definite  dependence 
of  on  2:,  .  Consequently,  given  a  2*  >  the  total  pressure  drop  across  a  bow 

wave  experienced  by  a  streamline  passing  through  the  shock  wave  at  2*  can  be  deter¬ 
mined.  Now  at  any  station  that  is  being  analyzed  the  static  pressure  at  the 
edge  of  the  boundary  layer  is  known.  From  the  total  and  static  pressure  at  the 
edge  of  the  boundary  layer  and  the  free-stream  total  enthalpy,  the  static  tempera¬ 
ture  at  the  edge  of  the  boundary  layer  can  be  determined.  As  yet  the  total  pres¬ 
sure  has  not  been  defined  as  a  unique  function  of  Xn  (i.e.,  the  relationship  be¬ 
tween  2*  and  X,,  is  not  yet  known).  Consequently,  at  each  station  Xu  for  which 
the  analysis  is  carried  out,  a  curve  of  2*  vs  "Ti  must  be  drawn.  Such  a  curve  is 
effectively  one  equation  in  the  two  unknowns  2*  and  “TV  .  With  a  knowledge  of  TV 
and  ^  reference  26  can  be  used  to  define  V\s  such  that  Mr  can  be  obtained  from 

U*  *  I2  .  (4-45) 

Equation  (4-45),  in  conjunction  with  the  curve  of  2*  vs  Tg  is  sufficient  to  per¬ 
mit  a  curve  of  vs  Mg  to  be  drawn  at  each  station.  Continuing  in  basically 
the  same  manner,  and  employing  reference  26,  it  is  possible  to  plot  a  curve  of 
vs  £,  and  vs  1,  at  each  station.  Such  curves  come  in  very  handy  when  a 
numerical  solution  is  being  obtained  by  use  of  a  desk  calculator. 

The  local  heat  transfer  Q  at  any  point  can  be  determined  from  the  following 
equation,  6“r 


(4-46) 
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Applying  the  Dorodnltzyn  transformation  and  converting  the  temperature  derivative 
to  an  enthalpy  derivative,  the  heat  transfer  at  the  vail  can  be  written 


(4-47) 


where  the  enthalpy  derivative  has  been  evaluated  using  the  enthalpy  distribution 
of  equation  (4-36). 


The  local  shear  stress  at  any  point  can  be  determined  by  considering  the 
equation 


(4-46) 


Evaluating  from  the  polynominal  expression  for  the  velocity  profile 

(equation  (4-10))  the  local  Bhear  stress  can  be  written 


(4-49) 


V.  SOLUTION  FOR  THE  ADIABATIC  WALL  CASE 

In  Section  IV  an  expression  for  Va  as  a  function  of  A  was  developed 
(equation  (4-43))  for  a  condition  of  heat  transfer  to  the  wall.  This  term  could 
also  be  developed  starting  with  the  total  enthalpy  distribution  applicable  for 
zero  heat  transfer  (equation  (4-33)).  The  integration  of  equation  (4-27)  leads 
to  the  following  expression  for  application  to  the  case  of  zero  heat  transfer: 


o».  Vu.-Vw ,&\-f 
y  — ^ i m) 


(5-D 


By  examining  the  energy  Integral  equation  it  can  be  shown  that  9*5  O  for  the 
adiabatic  wall  solution.  For  the  condition  of  zero  heat  transfer  the  term  involving 
the  integral  of  in  equation  (4-28)  is  identically  zero  and  hence  must 

also  be  zero,  ms  being  the  caBe,  equation  (3-1)  becomes  the  second  egression 
for  ,  the  first  being  equation  (4-39).  If  is  eliminated  from  both  equations 

the  following  equation  results, 


l-Y, 


u t 


Zy  -  a. 

4»o  TX* 


(5-2) 


where  ,  the  recovery  factor  at  the  wall,  has  been  defined  as 


^  cuts*  -  Vi* 


(5-5) 
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On  a  blunt  body,  such  as  the  one  being  analyzed,  the  pressure  gradient  parameter 
/V.is  always  positive.  This  is  shown  in  equation  (4-14)  whereA.  is  expressed 
as  a  function  of,  among  other  terms,  ■  The  term  fc  **/u)*,  is  positive 

everywhere  on  the  plate  (see  Figures  8  and  9)  as  are  the  other  terms  in  the  ex¬ 
pression  for  A  ,  consequently  A.  itself  must  always  be  positive.  A  plot  of  the 
recovery  factor  at  the  wall  is  shown  in  Figure  12  for  Prandtl  numbers  of  .67,  .72, 
and  .77. 

The  local  recovery  factor  \f  varies  through  the  boundary  layer  as  a  function 
of  where 

-  W 

V\vh*  (5J,) 

Substituting  the  enthalpy  distribution  for  zero  heat  transfer  into  equation  (5-4), 
the  following  expression  for  the  local  recovery  factor  is  obtained 


r  -  >Tur  ♦  O-ror)  K.to.)  ♦  ^a<-^  (5_5) 

where  equation  (5-1)  has  been  used  to  replace  fj,  in  equation  (4-53).  A  plot  of 
the  local  recovery  factor  for  .72  is  shown  in  Figure  13  for  .A.  of  -6,  0, 

and  6. 


Equation  (5-2)  can  be  considered  an  expression  between  ~TL. .  Tt  ,  and^V. 
Employing  it  in  the  system  of  equations  derived  in  Section  IV,  it  becomes  possible 
to  define  the  adiabatic  wall  temperatures  consistent  with  a  given  pressure  distri¬ 
bution.  Consequently,  the  integral  analysis  can  be  used  to  predict  the  surface 
temperature  distribution  for  flow  over  an  adiabatic  flat  plate  with  a  finite  pres¬ 
sure  gradient.  One  such  method  is  outlined  in  Appendix  VII. 


VI.  STAGNATION  POINT  ANALYSIS  FOR 
TWO-DIMENSIONAL  BLUNT  BODIES 


In  order  to  start  the  step-by-step  solution  of  the  system  of  equations  derived 
in  section  IV  it  is  necessary  to  evaluate  the  parameters  and  ©*  at  the 

stagnation  point.  The  momentum  integral  equation  for  a  compressible  fluid  can  be 
expressed  as 


f  -M*1i . 


(4-16) 


The  energy  integral  equation  as  developed  in  Appendix  I  (a  different  form,  but  a 
more  classical  expression  than  that  of  section  IV)  can  be  written 


i©I 

<AX 


JlrU* 


4 


(6-1)  * 
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Equations  (4-l8)  and  (6- l),  as  they  stand,  are  not  applicable  at  the  stagnation 
point  because  *4 g  ,  which  appears  in  the  denominator  of  both  equations,  is  zero  at 
a  stagnation  point.  Consequently  the  derivatives  take  the  indeterminate  form 


O 

~  5 


(6-2) 


The  onlv  possible  way  for  equation  (5-2)  to  have  a  finite  value  at  X  =  0  would  be 
for  and  <*€Ph/«JX  to  be  of  the  form  0/0.  The  values  A.  and  ti,,  which  force 

the  numerator  of  equations  (4-18)  and  ( 6— l)  to  go  to  zero  simultaneously  will  be 
the  stagnation  point  solution  for  A  and  .  With  the  values  of  A  and  so 

determined,  values  for  the  derivatives  in  equation  (6-2)  can  be  determined  in  the 
limit  (if  they  exist)  by  applying  L'Kospital's  rule. 


Two  equations  in  the  two  unknowns,  and  o ,  can  be  obtained  by  setting 
the  bracketed  terms  in  equations  (4- lo)  and  (6-1)  equal  to  zero.  Die  resulting 
equations  are 

a*]., 

(6-3) 

b...  =  A. (I")  , 

*  *  (6-4) 

z 

where  the  terms  have  been  dropped  because  they  are  zero  at  the  stagnation 

point. 


In  order  to  solve  equations^  (6-3)  and  (5-4)  for  and  -A-.an  analytic  ex¬ 
pression  must  be  developed  for  ST/&  .  An  expression  has  to  be  obtained  for  Tfr* 
in  terns  of  in  order  to  evaluate  the  integral  in  equation  (4-22).  This  choice 
cannot  be  an  arbitrary  one  since  an  expression  for”TYT*  has  already  been  defined. 
That  is,  we  actually  define  the  boundary  layer  temperature  profile  when  the  total 
enthalpy  (equation  (4-35))  and  velocity  distributions  (equation  (4-10))  are  defined. 
At  the  stagnation  point  where  V\4,»V»  ,  since  U*  is  everywhere  zero,  the  tempera¬ 
ture  distribution  is  only  a  function  of  the  total  enthalpy  distribution.  The  problem 
now  arises  of  finding  an  expression  forf/^Tji  as  a  function  of  ^h4  ,  so  that  the 
integral  of  equation  (4-22)  can  be  analytically  integrated.  Now  in  the  most  general 
case,  enthalpy  is  a  complicated  function  of  both  pressure  and  temperature  and  an  ex¬ 
plicit  relationship  T»TO\}  is  not  generally  available.  However,  it  is  always 
possible  to  arrive  at  some  good  approximate  expression  for  T-TCVt')  when  certain 
restrictions  are  placed  on  the  range  of  temperature  and  pressure  variation. 


An  approximate  expression  of  this  type  has  been  developed  in  Appendix  II.  Die 
resulting  equation  has  the  following  form 


X  =  .<1*4  +.o6 
Tt  h. 


(5-5) 
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and  has  the  following  restrictions: 


1. 

*  £  znocT  R  . 

2. 

3t*<? 

4  T  £■  2*700*  R  . 

^  P«- 

3. 

-ooi 

i  -p  £  lo 

<*k. 

(6-6) 

Applying  equation  (6-5)  to  the  integral  in  equation  (4-22)  and  then  substi¬ 
tuting  for  the  enthalpy  ratio  the  enthalpy  profile  function  (equation  (4-55)), 
there  results  the  following  expression  for  ,  valid  only  at  a  stagnation 

point, 


'  ur 


Substituting  for 
equation  (6-5)  becomes 


& 


a  from  equation  (6-?)  and  for  % 


(6-7) 

from  equation  (4-20^ 


.l%8 


rto 


1Z 


]. 


(6-8) 


Equation  (6-4)  for 


can  be  rewritten  as 

A. 

\AA 


H? 

4  do 


302^) 


where 

tuted 


^/a 


(6-9) 


^  (equation  (4-45))  evaluated  at  the  stagnation  point  has  been  substi- 
for  the  term  in  equation  (6-4).  It  might  be  mentioned  that  the 

simplest  form  of  the  equation  formed  by  eliminating  b.,#  from  equations  (6-8)  and 
(6-9)  will  be  a  seventh  order  polynominal  in  A.  ,  and  even  then  the  term  involving 
fcC^Twas  considered  small  compared  with  one  in  order  to  simplify  the  resulting 
expression. 


Equations  (6-8)  and  (6-9)  do  not  generate  the  final  value  for  A*  and  idii. 
The  final  values  are  obtained  by  an  iteration  process  as  follows.  First  solve 
equations  (6-3)  and  (6-9)  foriL*  and  b,,4  .  Using  these  values  for  A*  and  bH# 
calculate  from  equation  (4-22).  The  temperature  ratio  in  equation  (4-22)  is 

determined  from  the  enthalpy  distribution  of  equation  (4-42)  as  a  function  of 
and  A.  as  determined  from  equations  (6-8)  and  (6-9).  This  value  of  5)^  is 
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f. 


substituted  into  equation  (6-3)  which  then  becones  the  equation  from  which  m 
is  determined.  With  this  value  of  equation  (6-9)  can  be  used  to  solve  for 
feii0  .  lhe  iteration  process  can  be  continued  until  the  values  for  fea,  end 
it^re  obtained  to  the  accuracy  desired. 

With  and  obtained  such  that  is  in  the  form  0/0,  it  is  possi¬ 
ble,  by  the  method  of  L' Hospital,  to  evaluate  the  limit  if  one  does  exist.  If,  in 
equation  (4-l8),  F,IX)  is  defined  to  be 


F.tx) 


(6-10) 


then,  taking  the  derivative  of  the  numerator  and  denominator  of  equation  (4-l8) 
with  respect  to  X,  the  following  equation  is  obtained 


a  (?) 


(5-11) 


1  Q* 

The  above  expression  is  used  to  determine  the  initial  slope  •  First,  how¬ 
ever,  the  terms  and  Utf.  must  be  evaluated.  The  term  V,*£»0#  is  defined 

in  Appendix  III,  and  can  be  reduced  to  a  relatively  simple  form  for  a  given  set  of 
conditions.  The  term  can  be  written 


(  i-12) 


and  its  derivation  will  be  found  in  Appendix  IV. 


Also  of  interest  at  the  stagnation  point  is  the  transformed  boundary  layer 
thickness  /X  .  This  can  be  obtained  from  equation  (4-15)  where  the  term 
is  obtained  in  Appendix  IV  as 


u. 


m 


(-J0. 


(5-13) 


Substituting  equation  (5-13)  into  the  equation  (4-15),  the  expression  for  the  stagnation 
point  boundary  thickness  in  the  transformed  plane  is 


A  ~  7/ 

°  * - g  /  (5-lM 

Once  A  has  been  determined  at  the  stagnation  point,  the  momentum  thickness  0 
can  be  obtained  from  equation  (4-23).  The  heat  transfer  at  the  wall  as  derived 
in  Section  IV  is 


A 


(4-47) 
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and  is  applicable  at  the  stagnation  point  when  the  expression  for  is  used. 
Making  the  required  substitutions 


ft,  ^ 

CRTO'  J 


1.08 


•jLjl 


(6-15) 


As  a  check  of  the  stagnation-point  analysis  two  flight  conditions  and  two  tunnel 
conditions  were  analyzed.  Table  I  summarizes  the  results  of  the  analysis.  The 
analysis  of  Case  I  was  plotted  in  Figure  14  along  with  the  results  of  references 
27,  28,  and  29.  The  over-all  agreement  between  the  integral  analysis  and  the 
above  references  is  encouraging. 

The  parameter  b%/Aj  appearing  in  equation  (6-15)  is  basically  the  only 
unknown  quantity  in  that  equation;  the  other  terms  will  be  known  at  a  given  con¬ 
dition.  The  parameter  /A*  is  plotted  in  Figure  15  against  Wvh.,  .  All 
the  points  analyzed  fall  on  one  curve.  This  curve  can  be  used  in  predicting  -VPr?,^ 
without  resorting  to  the  solutions  of  equations  (6-3)  and  (6-9)  to  determine  ® 
and  A.  .  This  curve  represents  a  rapid  method  for  estimating  two-dimensional 
stagnation-point  heat  transfer. 


VII.  DISCUSSION 

The  analysis  just  presented  is  adequate  to  cover  a  wide  range  of  conditions. 
However,  there  are  some  conditions  to  which  it  cannot  be  applied  directly,  and 
others  to  which  it  cannot  be  applied  at  all.  The  conditions  at  which  the  present 
analysis  is  limited  in  scope  can  be  considered  to  stem  from  two  basic  problem  areas: 

1.  Presence  of  strong  streamwise  pressure  gradients. 

2.  Presence  of  vorticity  in  the  free- stream  flow. 

First,  consider  the  problem  of  the  strong  streamwise  pressure  gradients.  In 
defining  the  velocity  profile  by  the  fourth-order  polynominal  in  i\,  there  is  a 
built-in  limitation  on  A  .  For  values  of  A  greater  than  12  the  velocity  profile 
overshoots  or  exhibits  a  velocity  somewhere  in  the  boundary  layer  which  exceeds 
the  free-stream  value.  This  is  physically  unlikely,  even  under  the  strongest  favor¬ 
able  pressure  gradient.  Yet,  in  analyzing  high  Mach  number  conditions  there  are 
many  times  when  the  solution  to  the  system  of  equations  is  possible  with  12, 

but,  because  of  the  Inherent  limitations  on  A  ,  the  solution  to  the  problem  must 
stop.  Unfortunately  pressure  gradients  strong  enough  to  cause A  to  exceed  12  are 
possible  in  many  situations,  such  as  at  the  leading  edge  of  blunt  plates  in  hyper¬ 
sonic  flow,  and  consequently  occur  in  regions  of  significant  interest  to  warrant 
an  attempt  at  a  solution  within  the  framework  of  the  basic  integral  analysis. 
Fortunately  such  an  approach  does  exist.  The  details  of  this  solution  will  be 
found  in  Appendix  VI  along  with  a  summary  of  the  pertinent  equations.  Only  the 
basic  features  of  the  analysis  will  be  discussed  here. 

The  significant  step  in  obtaining  a  solution  consists  in  a  redefinition  of 
the  basic  velocity  profile.  Steiger  (32)  defines  a  velocity  profile,  which  satis¬ 
fies  all  the  boundary  conditions  of  equation  (4-6),  in  the  following  manner. 
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(7-1) 


When  this  profile  is  used  in  the  analysis!  there  are  no  limitations  of  the  magni¬ 
tude  ot  Cl  except  that  A  *2  .  Since  this  is  being  proposed  as  an  alternative 
solution  for  flow  fields  with  strong  streamvise  pressure  gradients,  this  is  ef¬ 
fectively  no  restriction  because  streamvise  gradients  sufficient  to  make 
have  values  of  Ou>A  .  The  fact  that  the  problem  of  strong  favorable  pressure 
gradients  on  cylindrical  leading  edges  is  more  severe  the  higher  the  ffach  number 
can  be  shovn  in  the  following  manner. 


The  data  of  Figure  6,  for  cylindrical  leading  edges,  illustrate  the  indepen¬ 
dence  of  the  pressure  distribution,  when  plotted  as  P/p.  ,  and  free-stream  Mach 
number.  This  being  the  case,  the  maximum  pressure  gradient  expressed  in  the  form 
must  also  be  independent  of  the  free-stream  Mach  number  and  the  term 


JyVfowhlch  is  used  in  the  analysis  as 


can  be  obtained  from 


d?  -  4  F> 
d<t>  ♦  * 


(7-2) 


Now  increases  with  Mach  number  (for  a  given  free-stream  static  pressure), 
consequently  the  maximum  pressure  gradient  is  a  unique  function  of 

Mach  number.  Thus  the  problems  associated  with  strong  favorable  pressure  gradients 
will  become  more  severe  at  the  higher  Mach  numbers. 


It  might  be  noted  that,  contrary  to  the  published  results  of  reference  (52), 
there  is  no  unique  relation  between^,  and  Cl  ,  but  there  are  at  least  three  re¬ 
lationships  which  differ  in  all  but  one  respect.  The  three  equations  between  <X 
and  .A  can  be  obtained  by  evaluating  the  following  expressions: 


using  both  the  polynominal  profile  (equation  (4-10))  and  the  exponential  profile 
(equation  (7-l)).  Assuming  that  there  is  an  unique  expression  between  <X  and ./V  , 
then  the  above  three  expressions  should  be  invariant  with  respect  to  the  choice  of 
velocity  function.  Evaluating  the  above  three  expressions  using  both  velocity  pro¬ 
files  and  equating  like  terms  yields  the  following  three  equations: 


A= 


3A  -n 

OL  I 


(7-3) 
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2. 

A  =  . 

(7-4) 

3. 

-A.  =  clCo--*')  . 

(7-5) 

The  one  condition  common  to  all  three  equations  is  the  fact  they  all  pass  through 
the  point  CL  =  k,  A.  =  12. 

This  is  a  very  important  fact  when  solving  a  numerical  problem,  as  can  be 
illustrated  in  an  example.  Consider  an  analysis  which,  because  of  the  strong 
streamwise  pressure  gradient,  has  been  started  with  the  system  of  equations  de¬ 
veloped  in  Appendix  VI.  As  the  solution  proceeds  point  by  point  around  the  body, 
the  value  of  A  will  reach  a  peak  and  then  decline.  Anytime  (X  becomes  less  than 
4,  the  analysis  can  be  carried  on  by  either  method.  This  change  from  one  system 
of  equations  to  another  can  be  made  because  of  the  step-by-step  nature  of  the 
solution.  The  only  fact  of  any  consequence  at  station  JCn+i  is  the  results  of 
the  analysis  at  Xn  ,  not  the  system  of  equations  used  to  obtain  the  results.  In 
the  discussion  of  the  example  it  has  been  implicity  assumed  that  the  results  ob¬ 
tained  by  either  method  of  analysis  are  essentially  the  same  inthe  region  ^4  • 

The  second  problem  area,  vorticity  existing  in  the  free-stream  flow,  has  been 
investigated  by  many  authors  ( 33>  34,  35).  The  basic  facet  of  the  problem  is  that 
with  vorticity  in  the  inviscid  flow  field  the  boundary  condition  at  the  outer  edge 
of  the  boundary  layer  involving  is  no  longer  valid.  The  effects  of  vorticity 

become  appreciable  when  the  velocity  gradient  at  the  outer  edge  of  the  boundary 

layer  is  of  the  same  order  of  magnitude  as  (32). 

& 

The  region  on  the  body  where  the  inviscid  velocity  gradient  becomes  important 
is_  that  point  on  the  plate  where  the  streamlines  entering  the  boundary  layer  have 
passed  through  the  highly  curved  portion  of  the  shock  wave.  At  low  Reynolds  numbers 
this  occurs  in  the  stagnation  point  region,  but  at  higher  Reynolds  numbers  this 
could  occur  on  the  aft  portion  of  a  long  body.  The  vorticity  problem  is  inherently 
more  severe  at  the  leading  edge  region.  Fortunately  at  the  leading  edge  of  a  two- 
dimensional  blunt  flat  plate  the  vorticity  interaction  becomes  a  problem  about  the 
same  time  that  the  assumption  ^  «  o  breaks  down.  Consequently  the  entire  analy¬ 
sis  is  invalid  and  there  is  little  hope  for  a  solution  within  the  framework  of  an 
integral  approach.  The  applicability  of  this  analysis  can  always  be  checked  by 
comparing  the  value  of  S  at  the  stagnation  point  with  the  leading  edge  radius  h  . 

If  Mr  «  |  ,  then  the  integral  approach  as  formulated  in  this  analysis  will  not 
apply. 

Now  the  other  region  where  vorticity  can  become  a  problem,  on  the  aft  ends  of 
long  plates,  can  be  treated  within  the  framework  of  the  integral  analysis  in  the 
following  manner.  For  a  constant  energy  flow  field  the  vorticity  at  the  outer  edge 
of  the  boundary  layer  can  be  expressed  as  a  function  of  the  entropy  gradient  normal 
to  the  streamlines  in  the  following  manner  (36), 


4- 
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This  co-ordinate  ki  is  normal  to  the  streamlines.  On  the  aft  end  of  very  long 
bodies  the  streamlines  are  very  nearly  parallel  to  the  body  surface  and  conse¬ 
quently  gradients  normal  to  the  streamlines  in  that  vicinity  would  not  differ 
greatly  from  gradients  normal  to  the  body.  In  terms  of  the  entropy  gradient 
normal  to  the  body  the  vorticity  at  the  outer  edge  can  be  expressed  as 


(7-7) 


For  two-dimensional  flow  the  vorticity  can  also  be  written 


(7-8) 


where  the  term  ,  therefore  can  be  expressed  as 


(7-9) 


In  the  analysis  presented  in  Section  IV,  was  illicitly  assumed  to  be 

zero  and  consequently  the  velocity  gradient  at  the  outer  edge  of  the  boundary 
layer  was  taken  to  be  zero.  With  vorticity  interaction,  or  entropy  gradients,  at 
the  outer  edge  of  the  boundary  layer,  equation  (7-9)  becomes  the  new  boundary  con¬ 
dition  on  the  velocity  profile.  The  second  derivative  will  remain 

zero  for  purposes  of  simplifying  the  analysis.  Hie  experimental  data  of  reference 
9,  would  seem  to  justify  this  assumption. 

Equation  (7-9)  defines  in  terms  of  the  entropy  gradient  at  the  outer 

edge  of  the  boundary  layer.  In  order  to  evaluate  this  derivative  an  entropy  pro¬ 
file  of  the  form 


X 


(7-10) 


is  assumed,  where  A  is  the  co-ordinate  normal  to  the  body  such  that, 


(7-11) 


The  subscript  &  refers  to  conditions  at  the  shock  wave.  The  boundary  conditions 
on  equation  (7-10 )  are: 


A-i 
A  - 1 
A=o 


A  -  1 

* 

- 

5  „  Sr 
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5. 


(7-12) 


The  entropy  behind  the  shock  wave  at  the  station  in  question  can  be  found  from 
the  thermodynamic  parameters  of  the  fluid  at  that  point.  The  derivative 
at  X»  \  can  be  determined  from  the  shock-wave  shape  and  the  free-stream  Mach* 
number.  A  good  representation  of  the  shock  wave  downstream  of  the  leading  edge 
region  is  given  by  the  blast-wave  theory  (l*v,  15,  16,  17,  18,  19)  and  there  is 
good  experimental  verification  of  its  usefulness  (37  ,  58)  •  The  third  boundary 
condition  is  connected  to  the  basic  boundary  layer  solution  in  that  it  is  a 
function  of  the  parameter  .  The  fourth  boundary  condition  is  developed 

from  a  conservation  of  entropy  analysis.  The  streamline  ,  in  addition  to 
having  a  constant  until  it  enters  the  boundary  layer,  is  also  isentropic. 
Every  other  streamline  crossing  the  shock  wave  with  will  cross  the  plane 

normal  to  the  body  at  ***  with  the  same  entropy  it  had  after  it  crossed  the  shock 
wave.  This  represents  a  certain  fixed  entropy  at  the  plane  normal  to  the  body  at 
X*  .  Equating  this  fixed  sum  %%  to  the  integral  of  the  assumed  profile  allows 
the  fourth  coefficient  to  be  determined  from 


£  g 


(7-13) 


The  entropy  behind  the  shock  wave,  i  & ,  and 
in  the  plane  normal  to  the  body  at  ,  can 
the  following  equation, 


^  S  -  S 


p  uiy 


consequently  the  entropy  of  the  flow 
be  determined  at  the  shock  wave  from 


Sc«p  <1  2 
4 


( 7-1*0 


The  terms  ,  and  are  functions  of  flow  conditions  on  the  down¬ 

stream  side  of  the  shock  wave  which  can  be  determined  once  the  free-stream  con¬ 
ditions  and  shock-wave  shape  are  known.  The  integration  indicated  by  equation 
(7-1*0  can  be  carried  out  numerically  or  analytically  depending  on  whether  the 
ratio  of  specific  heats  is  the  same  on  both  sides  of  the  shock  wave.  Consequently 
it  is  possible  to  determine  the  coefficients  A,  of  equation  (7-10)  by  equating 
equation  (7-1*0  to  (7-13)  and  solving  for  Q-i  .  The  expression  for  the  velocity 
gradient  at  the  outer  edge  of  the  boundary  layer  when  the  vorticity  is  present 
(equation  (7-9))  becomes 


(7-15) 


Equation  (7-15)  can  then  be  used  to  define  It  might  be  noted  that  to 

obtain  a  numerical  value  for  ( for  any  given  problem  conditions  existing 
at  the  outer  edge  of  the  boundary  layer  at  must  be  known.  This  information 
would  not  be  generally  known  and  a  straightforward  solution  would  require  adding 
the  pertinent  equations  of  this  section  to  the  general  analysis  and  determining 
all  the  parameters  by  a  simultaneous  solution.  This  would  be  a  lengthy  procedure 
and  as  much  could  be  obtained  by  solving  the  zero  vorticity  problem  first  and  then 


25 


using  the  zero  vorticity  solution  as  the  starting  point  to  evaluate  the  effects 
of  vorticity  using  the  analysis  Just  presented. 


The  question  might  next  be  asked:  when  should  the  vorticity  correction  be 
applied?  The  answer  can  only  come  from  a  numerical  solution  of  a  particular  case. 
When,  in  the  course  of  a  solution,  the  Btreamlines  entering  the  boundary  layer  are 
crossing  the  shock  wave  in  the  highly  curved  portion  of  the  shock  wave,  it  .would 
be  well  to  examine  the  term  to  determine  if  it  is  of  the  order 

If  they  are  of  the  same  order,  the  vorticity  interaction  analysis  should  be  applied 
to  ascertain  what  effect,  if  any,  the  vorticity  at  the  outer  edge  of  the  boundary 
layer  has  on  boundary  layer  parameters. 

The  analysis  presented  in  this  paper  is  also  applicable  to  blunt  flat  plates 
at  angle -of -attack  greater  than  zero.  The  only  requirement  is  that  the  surface 
temperature  and  pressure  distributions  be  available.  At  non-zero  angle-of -attack, 
analytical  methods  for  predicting  pressures  over  the  body  are  very  unreliable  and 
it  would  be  expected  that  the  pressure  distribution  would  have  to  be  obtained  ex¬ 
perimentally.  An  analysis,  in  basically  the  same  form,  could  be  developed  to  de¬ 
fine  the  compressible  boundary  layer  characteristics  on  an  axi-symmetric  blunt 
body.  The  essential  features  of  such  an  analysis  would  be  the  same,  except  that 
for  an  axi-symmetric  body  the  analysis  would  not  be  applicable  at  angles -of -attack 
other  than  zero  because  such  a  flow  field  is  no  longer  two  dimensional. 


VIII.  CONCLUSIONS 

The  Von  Karman  momentum  integral  has  been  extended  to  include  prescribed  wall 
temperature  and  pressure  gradients.  The  effects  of  entropy  gradients  in  the  in- 
viscid  flow  have  been  considered  and  a  first-order  correction  to  account  for  en¬ 
tropy  gradients  has  been  developed. 

The  analysis  has  been  carried  out  for  blunt  flat  plates  with  cylindrical 
leading  edges  at  zero  angle-of-attack.  The  analysis  itself,  however,  is  valid 
for  any  reasonable  shape  or  angle-of-attack  provided  the  pressure  distribution 
and  surface  temperature  distribution  are  known.  The  basic  features  could  be 
applied  to  develop  a  similar  system  of  equations  valid  for  axi-symmetric  bodies 
at  zero  angle-of-attack. 

The  problem  of  very  high  favorable  pressure  gradients  was  considered  and  an 
alternative  solution  developed  based  on  a  velocity  profile  which  does  not  exhibit 

profiles  in  which  U  >M.  . 

© 
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APPENDIX  I 


Derivation  of  the  Energy  Integral  Equation  for 
Use  at  the  Stagnation  Point 


If  equation  (4-28)  is  differentiated  with  respect  to  x, the  following  equation, 
can  be  written; 


(I-D 


Substituting  for  (-qj  from  equation  (4-47),  equation  (I-l)  may  be  written 


Ik  t, 

■»*  tV-A 

Equation  (l-2)  can  also  be  written 


(1-2) 


J>U 

J  "dX 


+u.|f)  = 


The  compressible  Bernoulli  equation  may  be  expressed  as 


(1-3) 


U 


da  +  if 

T 


*  o 


and  from  the  definition  of  the  speed  of  sound  it  is  possible  to  write 


Combining  the  last 


two  equations  in  such  a  way  as 

if  .  du 


to  eliminate 


there  results 


(1-4) 


After  substituting  equation  (l— 4)  into  equation  (1-3)  and  rearranging  the  terms, 
one  obtains  the  following  equation 


d©. 

A  X 


TCjkS  A 


A 


(1-5) 
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Then;  if  both  sides  of  equation  (1-5)  are 
equation  is  obtained 


aeC-  _  T*» 


multiplied  by  ,  the  following 


\l*  l  %  ^ 

-£)  '  ^  has  replaced  ,  with  A  determined  from  equation  (4-15) 
and  the  Bernoulli  equation.  Factoring  2  from  both  terms  on 
the  right  hand  side  of  equation  ( X— 6) ,  one  obtains  the  equation 


2A*£T>* 

ax’" 


*  l*sr* 


♦ 


(1-7) 


Equation  ( 1—7 )  is  used  in  conjunction  with  equation  (4-l6)  to  define  the  stagnation 
values  for  and  Ji.  . 
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APPENDIX  II 


Approximate  Relationship  Between  — 

Tjj 


and 


In  evaluating  the  terms  A  and  at  the  stagnation  point  the  term  % 
appears  in  the  analysis  in  such  a  way  that  it  is  necessary  to  evaluate  the  integral 
•  The  assumed  enthalpy  profile  implicitly  defines  the  ratio  X  hut 
an  explicit  expression  between  and  is  not  available,  consequently  an 

attempt  is  made  to  obtain  an  approximate  relationship. 

The  expression  for  the  relationship  betveenX^Tx'  and  was  obtained  from 

a  correlation  of  the  data  appearing  in  reference  (26)  and  reference  (31).  The  data 
were  analyzed  in  the  following  manner.  First,  for  assumed  values  of  “T>  equal  to 
3oO°R,  900°R,  ldOO°R,  and  2700°R,  curves  ofT/-^  Vs  V)  ^  were  drawn  for  pressures 
of  .001  *  p'e.,*  10.  It  was  observed  that  for  each”l^  there  was  a  region  of  the 

curve  which  was  independent  of  the  pressure  and  which  could  be  adequately  represented 
by  a  straight  line,  of  the  form  x‘rf'AX*«  • 

The  values  of  the  constants  A  and  B  which  best  fit  the  data  of  the  real  gas 
tables  are  A  =  .92,  and  E  =  .0o.  The  curve 


T 

T* 


(ii-D 


along  with  data  from  references  26  and  31  are  plotted  in  Figure  lo.  As  can  be 
seen  from  Figure  10,  the  curve  represented  by  equation  (il-l)  is  a  valid  approxi¬ 
mation  for  a  relationship  between and  M,*  •  The  expression  remains  valia 
as  long  as: 


lb/ft2. 


1. 

.  Oo  |  6r 

-3* 

a*. 

»o 

2. 

2>foo  < 

T> 

£ 

2  7oO*£ 

3- 

T 

< 

2noo°f? 

( H-2 ) 

surface 

pressure  is 

referenced  to  sea  level 

standard  pressure  of  2116 
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APPENDIX  III 


Let 


where 


f6V.\ 

Derivation  of  1 ——  1 

* 

be  defined  as 

Ffcx). 

(m-i) 

1. 

2. 

fi,  cx)  -  Tm-' 

3. 

6.  (**>.)*  -| 

(HI-2) 

The  term  K  is  a  new  variable  defined  as, 

*  - 


or 


4.  G^OO  *  -/i.  • 


(ni-3) 

(IIX-4) 


Using  the  definition  of  fx*)0  (equation  (lll-l)),  its  derivative  with  respect  to 
X, applicable  only  at  the  stagnation  point,  can  be  written 


where  it  has  been  implied  in  equation  (III-5)  that: 


1.  K  «•  <jaCX>  • 

2.  b„  .  (ui-6) 
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Each  term  in  the  expression 


1.  (&) 
(,-bXJc 


[I*]. 


will  be  investigated  separately. 


(g) * 


(in-?) 


b&t  =r  -  ^ 

ax  T5  dX 


.  (§\'  /3%) 

ax  Ux  J 


(m-8) 


Also  from  equation  (6-7) 


LA].  L  **©  X,  0  7s~  V  aT  gT  - 


(III-9) 


from  which 


[tf];  «&].*  *  v 

Substituting  equations  (III-8)  and  (III-10)  into  equation  (III-7)  one  obtains 

The  term  l^x].  will  come  from  the  known  stagnation  point  wall  temperature 


The  term 
distribution 


will  come  from  the  known  stagnation  point  wall  temperature 


ISE1 

Lax  1 


(in-12) 


Iftl.  •  (ft),  fit). 


(HI-13) 
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and 


Now 


®,i±l 


-Aifr 


(III-14) 


and 


^Vdt*^eUHL. 


i  A 

Ti.e  term  5;  -  is  evaluated  from  equation  ( b-20 )  as, 

cix 


(m-15) 


I£%1  = 

ImI 

37-  A„  SAi 
3  l-M 

(in-16) 


Substituting  equations  (III-14),  (III-15)  and  (III-16)  into  equation  (III-13)  the 
expression  for  C*Vak>.  becomes 

.1  as  _  3  ^  5  ^ 

33  2  *“-48  0 


( 


=  C%\ 

3K  -'o  4> 


The  term 


B7  -  Aa  -  j?S  A* 

(  j0  in  equation  (III-12)  can  be  written 


(III-17 ) 


V®.  I(ff  -®  a], . 


(III-16) 


fUA. 

I  3*  ;0 

The  partial  derivatives  of  ^  and  ^  with  respect  to  <A  are  obtained  from 
equations  (III-9)  and  (III-16),  respectively,  and  comes  from  equation 

(HI-15).  When  the  proper  substitutions  are  made,  equation  (lll-l8),  for/d^\ 
can  be  written  l 


(»).■ 


5‘Ci-Sr) 


S7-  A  -  2?  A1 


(m-19) 


Substituting  equations  (III-17)  and  (III-19)  back  into  equation  (III-12)  generates 
the  following  equation  for  C«vfcO. 

(t£V  ^ 


1  ^ 
4 


*7  -  A  -  -S^A* 

■* 


1 


(III-PO) 


The  derivative  ja^/dV>„].  can  be  evaluated  by  differentiating  Equation  (III-9) 
with  respect  to  and  is  found  to  be 

(bft  ^  _  .138  t  _±6  W 

1 )~  •WrUcfm  . 

the  derivative  ^ 


With  the  above  value  for  the  derivative 
becomes 


4. 


dJCl 

ixl 


(III-21) 


Starting  with  equation  (4-23)  and  rearranging  it  such  that 


e*  JV  i«l, 

k*  a  % 


enables  the  derivative 


9t\ 


to  be  obtained  as 


1  *  if  T*1 

r  ,Kl  „  r  i  4?\  A  m- .  A-  ijgj  •) 
[£]  *  K-\_e*  Xx  as*  Ux)  /fat?  J. . 


The  term 


zero 


(I 11-22) 


The  term  V dXJo  is  the  pentameter  that  is  being  sought  and  the  term  4'4/dx 
is  identically  zero  at  the  stagnation  point. 


(111-23) 

x 


deserves  a  closer  examination. 
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First 


R  T„ 


(m-24) 


and. 

jX  -  c  Tu» 

S*  u r  -  ■==, 

•8 


yUf 

In  the  last  expression  the  terms 
only  convenient  reference  quantities, 
the  term  _ can  be  written 


P„tJ  .  VHr 
>7  '  O 


^  are  treated  as  constant, 
Using  the  above  definitions 


(III-25) 


since 

for  A  and 


are 

JU*, 


J.  /  4  _L 

'  d  X  )»  Tu,  <&  X 


The  term  tr.Ux  is  identically  zero  at  a  stagnation  point  and  consequently 
equation  (HI-25)  reduces  to 


(III-26) 


5. 


It  has  been  implicity  assumed  that  a  relationship  exist r  such  that 
This  relationship  is 

Vx,,  ' 


ftk.AT.r.  sS; 

i  26ZO 


(6-9) 


This  expression  results  from  setting  the  numerator  of  the  energy  integral  equation 
equal  to  zero  and  substituting  for  the  stagnation  point  value  of  .  Differ¬ 

entiating  equation  (6-9)  once  with  respect  to  X  one  obtains 


Now  the  term 


dfe,.  £«  iT  ♦  S  £ 
5^-  aT  a*  av 


(111-27) 


where  again  the  term  av 
can  be  written  as 


<»P„o 


at  the  stagnation  point.  Hence  the  term 


x'o 


(.a*V.  ^>T  <*x 

/Afc>„  \ 

After  the  terms  in  equation  (III-27)  are  rearranged  the  expression  fo  *  J , 

can  finally  be  written  as 


fdAr 

-L.  •*  / 

,  0*  J.. 

sk  Y1 

UxJ0 

i  <jx  1 

.l4A  vw1 

.  4  A 

(-  zf). vdc<’“(  8T 


A. 

144  • 

(III-26) 

/AAn 

Using  equations  (III-15)  and  (III-26),  Ux).  =an  be  found  as 

m-  . 

Combining  the  results  from  sections  1-5  of  this  Appendix  there  results  an 
expression  for  ^ which  has  the  form 

~  (III-30 ) 

where is  some  parameter  dependent  on  and  X  . 
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APPENDIX  IV 


Derivation  of 


The  compressible  Bernoulli  equation  can  be  written 


u4  4.  . .  i  4j> 

dx  J*s  AX  • 


(IV-1) 


The  pressure  distribution  applicable  to  a  blunt  body  as  obtained  from  the  modified 
Newtonian  theory  is 


if  -  cos  2X„  r  4'-  '7 

p,  L  p>  J  • 


( IV-2 ) 


Differentiating  equation  (IV-2)  with  respect  to  X  ,  the  expression  for  the  pres¬ 
sure  gradient  is  obtained  as 


(IV-5) 


Combining  equations  (IV-1)  and  (lV-3),  one  obtains 

^  r  &*.  ji 

fs  U,  J  ■ 


im.- 15 


( IV-4  ) 


Now  at  the  stagnation  point  both  sin4X,  and  go  to  zero  and  the  expression 
for  t*v**].  must  come  from  L'Hospital's  method  for  evaluating  indeterminate 
forms.  Using  the  method  of  L'Hospital,  the  expression  for  at  the  stag¬ 

nation  point  becomes 


At  the  stagnation  point  the  term 
derivative  to  be  written  as 


llH*  \ 


■4h]  ,  yrstc 

**L  ~2 


The  term 


KrJ 


(lV-i>) 

which  allows  the  velocity 

(IV-6) 


comes  up  in  obtaining  an  expression  for  A*  .  It 


can  be  evaluated  from  equations  (IV-1)  and  (IV-6)  as  follows: 


56 


(IV-7) 
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APPENDIX  V 


Method  of  Numerical  Solution 


The  equations  developed  In  Section  IV  can  be  solved  in  the  following  manner. 
First  equations  (6-8)  and  (6-9)  are  solved  by  trial  and  error  to  yield  the  stag¬ 
nation  point  values  for  b„  and  A.  .  Knowledge  of  b„  and  A.  enables  one  to 
calculate  A*  and  from  equations  (6-14)  and  (6-15),  respectively.  With  A* 

known,  it  is  possible  to  define  the  stagnation  point  non-dimensional  velocity  pro¬ 
file.  This  can  be  considered  as  a  limiting  profile,  since  at  the  stagnation  point 
i4s*o  .  With  b„  known,  the  total  enthalpy  profile  at  the  stagnation  point  can 
be  defined  in  terms  of  the  non-dimensional  transformed  co-ordinate  »\  .  Since  the 
pressure  at  the  stagnation  point  is  known,  the  temperature  profile  can  be  obtained 
from  real  gas  tables  (26)  and  the  enthalpy  profile.  Integrating  equation  (4-22) 
numerically,  using  the  temperature  profile  obtained  from  the  enthalpy  profile,  and 
substituting  for  the  value  found  from  the  solution  of  equations  (6-8)  and 

(6-9),  the  term  can  be  found. 

The  next  step  is  to  find  the  stagnation  point  values  for  ©  and  *°'/4X  . 

The  term  obtained  from  equation  (4-23)  with  A*  defined  in  equation  (6-14). 

The  term/Ts  obtained  from  equation  (6-ll)  where  the  terms  Pi4*).  and  can 

be  obtained,  once  A,  and  b„#  have  been  defined  for  a  given  set  of  conditions, 
from  Appendix  III  and  IV,  respectively. 

Once  the  stagnation  point  values  for  O  and  are  known  the  solution 

can  be  extended  to  the  next  step.  The  step-by-step  method  of  solution  is  necessary 
because  the  differential  equation  in  the  system  (the  momentum  integral  equation) 
does  not  have  a  solution  in  closed  form.  As  in  any  numerical  solution  to  a  dif¬ 
ferential  equation,  the  accuracy  of  the  final  solution  is  an  important  function  of 
the  step  size. 

The  value  of  ©  at  the  first  step  is  found  from  the  stagnation  point  solution 
using  the  following  equation. 


(V-l) 

For  the  first  step  the  subscript  "0-1  "  would  correspond  to  the  stagnation  point 
values.  Consequently  the  necessity  of  the  stagnation  point  solution  is  apparent. 
The  stagnation  point  solution  is  the  first  link  in  a  sequence  of  events  that 
eventually  leads  to  the  solution  for  the  boundary  layer  characteristics  on  a 
blunt  flat  plate.  At  the  station  at  which  the  analysis  is  being  carried  out 
the  parameters  .  T-r,  and  Pr  are  known.  With  these  parameters  essentially 

fixed,  C*  can  be  obtained  by  a  method  completely  different  from  equation  (V-l). 

The  different  method  begins  by  arbitrarily  choosing  a  value  of  .  This  estab¬ 
lishes  a  trial  value  for  V*#  ,  the  downstream  total  pressures,  according  to  the 
shock-wave  analysis  of  section  IV.  Assuming  ,  the  ratio  can 

be  computed  and  this,  along  with  the  total  temperature,  establishes  the  static 
temperature  consistent  with  the  trial  2*  .  The  static  temperature,  Tf  ,  is 
obtained  from  the  following  equation. 
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(V-2) 

As  it  is  written,  it  is  not  an  explicit  expression  for”TV  because  "Xf  is  buried 
in  the  exponent  But  "Ti  can  always  be  obtained  from  (V-2)  by  a  trial 

and  error  approach.  Once  a  trial  "T*  has  been  established,  the  static  enthalpy 
and  hence  the  velocity  can  be  obtained.  The  static  enthalpy  is  obtained  from  a 
set  of  real  gas  tables  consistent  with  the  trial  “Tf  and  the  which  is  fixed 
by  the  choice  of  the  station  to  be  analyzed.  The  trial  velocity,  consistent  with 
the  initial  choice  of  and  arbitrary  choice  of  2*  ,  is  obtained  from  equation 
(4-45). 

At  this  point  the  equation  involving  the  continuity  of  mass,  equation  (4-ll), 
is  introduced  into  the  analysis.  The  expression  for  A*  ,  equation  (4-15),  is 
substituted  into  equation  (4-ll),  and  the  resulting  equation  becomes  an  expression 
for./\.  as  a  function  of  the  initial  choice  of  3.  ,  and  the  trial  values  of  Tr  and 

consistent  with  this  choice  of  2*  .  A  trial -and-error  solution  of  the  re¬ 
sulting  equation  will  yield  a  value  for /I  consistent  with  the  initial  choice  of 
Zm  .  With  -A.  so  determined,  ©can  be  determined  again,  this  time  from  equation 
(4-25).  This  value  of  is  compared  with  the  value  obtained  from  equation  (V-l). 
Different  values  of  are  chosen  and  the  whole^ process  repeated  until  a  value 
of  2,  is  chosen  such  that  the  trial  value  of  &  matches  that  obtained  from  equa¬ 
tion  (V-l).  At  this  point  the  value  of  2*  consistent  with  ,  such  that  the 
mass  flow  between  Vo  and  Hr.  is  equal  to  the  mass  flow  in  the  boundary  layer  at 
Xn  ,  has  been  determined.  Also  the  temperature,  total  pressure,  ratio  of  specific 
heats  and  velocity  at  the  outer  edge  of  the  boundary  layer  has  been  established. 

Next  the  value  for  is  determined.  This  is  done  by  numerically  inte¬ 

grating,  by  some  approximate  technique,  the  heat  transfer  between  the  stagnation 
point  and  .  Evaluation  of  dX  allows  °*/&  to  be  evaluated  from 

equation  (4-28).  With  this  value  of  »n4  equation  (4-1+5)  becomes  an  expression 
from  which  the  parameter  can  be  obtained.  With  this  value  of  and  the 
value  of  -A.  already  obtained,  the  heat  transfer  at  the  wall  can  be  obtained  from 
equation  (4-1+7). 

The  term is  next  used  to  define  the  velocity  profile  using  equation  (4-10), 
and  both  A  and  b,,  are  used  to  define  the  total  enthalpy  profile  from  equation 
(4-42).  Both  profiles  will  be  in  terms  of  the  transformed  co-ordinate  r\  .  To 
obtain  the  profiles  in  the  X«^  plane  the  static  temperature  distribution  is  required. 

The  static  temperature  distribution  is  obtained  from  the  total  enthalpy  and 
velocity  distributions  using  the  energy  equation, 

h  *  Vm  . 

(v-5) 

With  the  static  enthalpy  determined  from  equation  (V-5),  the  static  temperature  can 
be  obtained  from  a  set  of  real  gas  tables  as  a  function  of  and  The  ex¬ 

pression  between  U  and  is  obtained  from 
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(V-4) 


Equation  (V-4)  can  be  numerically  evaluated  to  obtain  a  tabular  relationship 
between  ^  and  .  This  technique  will  have  to  be  repeated  at  each  station 
for  which  plots  of  velocity,  or  temperature  in  the  X_<£  plane  ai*£  required. 

The  value  of  at  «\  •  >  i by  definition  £  and  is  exactly  the  integral  appear¬ 
ing  in  the  expression  for  sv*  .  Hence,  at  each  station  at  least  the  integral 
f1  T  must  be  evaluated  in  order  to  evaluate  SVa.  from  equation  (4-22). 

Tt  1  , 

With  S/A  evaluated,  the  expression  for  the  derivative  »7jx  is  now  com¬ 
puted  for  the  station  in  question  from  equation  (4-l8).  The  term  in  equa¬ 
tion  (4-l8)  is  evaluated  from  the  values  of  Uj  ,  ,  and.  ~Tg  obtained  in  the 

beginning  of  the  analysis. 


At  this  point  the  boundary  layer  characteristics  at  the  first  station  are 
known.  The  entire  analysis  is  then  repeated  at  the  next  station  starting  with 
equation  (V-l)  where  this  time  the  subscript  n-1  would  refer  to  the  values  of 
0'  and  4®Vix  just  obtained.  The  analysis  can  be  continued  until  the  entire 
body  has  been  analyzed. 


APPENDIX  VI 


Derivation  of  the  Equations  for  Use  with  Strong  Favorable 
Streamvise  Pressure  Gradients 


Start  by  defining  a  new  velocity  function  such  that 


_U  _ 


i  -  O-af. 


[4-10]*  (VI-1) 


Equation  (VI-l)  satisfies  all  the  boundary  conditions  of  equation  (4-8)  and,  in 
addition,  -©  where  i»  2  3  .  This  is  a  much  stronger  restriction  on  the 

profile  at  the  ouner  edge  of  the  boundary  layer  than  previously  used.  Evaluating 
the  integral 


\ 


I 

© 


0. 

Q.+  » 


(VI-2) 


which  is  used  in  equation  (4-6),  the  equation  for  2*  f rom  the  conservation  of  mass 
principle  becomes 


J>,U,  4  * 

*  *  <x*'  ■ 


[4-n]  (VI-3) 


x 

Utilizing  equation  (VI-l)  in  defining  A  from  the  momentum  equation  evaluated  at 
the  wall  (equation  (4-9))  the  new  expression  for  A*  results: 


A*=  ex  0-0 

^ 

AX 


[4-15]  (VI-4) 


Introducing  the  velocity  profile  into  the  compressible  momentum  equation 
(equation  (4-16)),  the  following  expression  results: 
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[4-18] 


(VI-7) 


At  this  point  it  might  be  stated  that  there  are  three  relationships  between 
A  and  Ol  .  The  result  is  that,  whenever  A  appears  from  C****V^Hl )  w  ,  substitute 
Q.0O.-I)  >  whenever  A  appears  in  the  form  from/«^Mr  >  substitute 

Ot  >  and  whenever  the  term  |  t  1  appears  from  Sl%ii^Arv.  >  substitute 

.  These  are  valid  substitutions  and  adhering  to  them  will  convert  all 
the  equations  in  the  main  body  of  this  report  to  this  new  velocity  profile  of 
equation  (VI-l). 


♦Numbers  in  square  brackets  refer  to  the  corresponding  equation  in  main  body  of  the 
analysis. 


/Ml 

Having  already  evaluated  J.  using  equation  (VI-l),  the  expression 

for  & */&  becomes 


[4-22]  (VI-8) 


Applying  the  Dorodnitzyn  transformation  to  the  definition  of  ^  (equation  (4-19)) 
and  then  evaluating  the  integral  using  the  new  velocity  profile  of  equation  (VI-l), 
one  obtains  the  following  result: 


&  -  - 

4 


[4-20]  (VI-9) 


The  parameter  Sz  appearing  in  the  profile  for  the  adiabatic  wall  total 
enthalpy  distribution,  can  simply  be  converted  by  recognizing  that  \  »  at  or, 


vh,-  w 

IfU Vtctov 


[4-59]  (VI-10) 


The  term  does  not  change  and  is  essentially  zero  when  **  ® 

The  biggest  change  in  the  whole  analysis  occurs  in  evaluating  .  It  is 

a  tedious  but  straightforward  procedure  and  when  completed  the  following  expression 
for  ,  applicable  to  heat  transfer  conditions,  results : 

q  h*,- W  \  3  _  +  _©  _  _3  -l 

'jf  ~  s  "a43  a+4 


- 


ZO  Oi\  a-*-4 


r,  p  1  **  ^4,-  W  \  i  I  .  ?  i  1 
t-  ***  h*,  bo  a-**  a+4  <*-♦*  j  . 


1 

cx+s  J 


[4-45]  (VI-11) 

The  recovery  factor  for  an  adiabatic  wall  can  be  obtained  from  equation  (VI-ll)  by 
equating  b|(  »  o  and  -  O  as  was  done  previously.  The  resulting  expression  is 

-L  _  JL  ♦  jg.  .11 

,  r  n  Off  A  0*5 

I-**,  -  a*  L1-  M - - - 

3  -  JSL  1  £-  .  3- 

5  0+4  Qf{ 

[5-2]  (VI-12) 
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The  stagnation-point  analysis  as  previously  presented  should  be  adequate  for 
all  pressure  gradients  since  at  the  stagnation  point.  For  comparison, 
the  equations  to  be  solved  for  and  will  be  given  below, 


Q.*4  cr„«-C  J 


£  AoCQ.-Ol  i  +  -  -3- 

c  J  20  Qo* » 


[6-9]  (VI-13) 

and 


\  \  \ 

m  V\+( 

/££?) 

3 

1^,  .1  I 

L 

V  'bTJ^ 

Oo*  l 

(VI-14) 

The  equations  presented  above  are  proposed  for  use  when  the  basic  analysis 
indicates  that  a  solution  lies  in  the  regime  where  .  This  analysis  of  the 

strong  gradient  case  using  the  equations  presented  in  this  Appendix  will  never 
generate  velocity  profiles  where  Cl  » ikg  .  It  is  limited  to  conditions  where  the 
profile  parameter  CL  is  greater  than  2. 


APPENDIX  VII 


Method  for  Predicting  Adiabatic  Wall  Temperatures 


The  adiabatic  wall  temperature  can  be  determined  within  the  framework  of  the 
Von  Karman  integral  analysis.  The  following  three  conditions  can  be  considered 
as  applicable  for  the  flow  over  an  adiabatic  wall: 


1. 


3. 


O.X  *  °  • 

[!•} 


»  O 


For  adiabatic  flow  over  a  plate  V\*»  *  V*#*,  and,  since  — yv>|  ,  it  is 

obvious  that  must  be  identically  zero.  In  addition,  if  there  is  no  heat 
transfer  to  the  wall  at  any  point,  then  the  heat  transfer  integrated  over  the  plate 
must  also  be  zero.  The  condition  follows  directly  from  the  second 

condition  and  equation  (4-20).  The  requirement  that  there  be  no  heat  transfer  to 
the  wall  eliminates  the  energy  equation  from  the  system  of  equations  to  be  solved. 
In  its  place  the  expression  ^  *  (equation  (5-2))  is  used. 


The  adiabatic  wall  temperature  can  be  defined,  consistent  with  a  prescribed 
pressure  gradient,  in  the  following  manner. 


1.  Solve  equation  (6-3)  for  (noting  ^»,0S  O  ). 

2.  Assuming  that  ,  use  equation  (6-ll)  to 

determine 

3.  with  A  •  determined  from  equation  (6-3)  and  the  (%). 

used  in  calculating  ,  evaluate  from 

equation  (4-23). 

4.  Pick  a  reasonable  value  for  tW ]  n  and  complete  the 

analysis  of  Appendix  V  to  determine  -/i  and  U*  consis¬ 
tent  with  the  choice  of  . 

5.  Using  the  value  of  .A  ,  and  U*  obtained  in  step  4, 

calculate  from  equation  (5-2). 

6.  Compare  the  value  of  obtained  in  step  5  with 

the  value  of  assumed  in  step  4. 
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7. 

8. 


Repeat  steps  4  and  5  at  the  same  station  until  the  value 

obtained  in  step  5  is  the  same  value  as  assumed 

in  step  4. 


Calculate 

(4-18). 


at  the  first  station  from  equation 


9. 


Find  Q  at  the  next 


step  from  equation  (V-l). 


10.  Start  back  at  step  4  and  repeat  the  analysis  until  ~Hw 
and  at  this  step  and  Ox  at  the  next  step  are 

known. 


In  the  above  manner  the  adiabatic  wall  temperature  may  be  determined  on  any  body 
over  which  the  integral  analysis  may  be  applied. 
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Table  1.  Summary  of  Stagnation  Point 
Heat  Transfer  Calculations 


CASE  I 

Altitude1=200,000  ft  Velocity* 10, 000  ft/ sec, r  =  12  in 


Vhti 

.429 

.215 

.107 

.054 

«*0 

5.77 

9.07 

11.54 

15.40 

CASE  II 

Altitude*  250,000  ft 

.215 

Velocity=15,000  ft/sec,r=12  in 

.107  .05k 

.027 

*-0 

15.70 

19.56 

25.68 

25.42 

CASE  III 

515  psia 

2060°R 

M1=l 2.26 

r  =  0.5  in 

Vhtl 

.659 

.621 

•  971 

1.000 

^0 

1.15 

.56 

.07 

0.00 

CASE  IV 

715  psia 

1960°R 

(■^=12.84 

r  =  0.575  in 

hw/hti 

.244 

- 

- 

- 

4v0 

10.7  k 

- 

- 

- 

^Atmosphere  and  thermodynamic  data  from  Reference  50. 


FIGURE  I.  PRESSURE  DISTRIBUTION  OVER  A  CYLINDRICAL 
BODY  AT  M,  ■  6.SS 
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FI8URE  2.  PRESSURE  OISTRI8UTION  OVER  A  CYLINDRICAL 
BODY  AT  M,- 12.28 
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FIGURE  3.  EFFECT  OF  ANGLE  OF  ATTACK  ON  THE 
PRESSURE  DISTRIBUTION  OVER  A 
CYLINDRICAL  LEADING  EDGE 
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FIGURE  4.  HYPERSONIC  BLAST  WAVE  ANALOGY 
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FIGURE  6.  CORRELATION  OF  THE  PRESSURE  DISTRIBUTION 
OVER  CYLINDRICAL  LEADING  EDGES 
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FIGURE  8.  PRgSS^JjlE MG,R,^|NT  0VER  A  CYLINDRICAL 
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FIGURE  9.  PRESSURE  GRADIENT  OVER  A  BLUNT  FLAT  PLATE  AT  M,*I2.2 


FIGURE  10.  SHOCK  WAVE  PROFILES  AT  M,« 12.28 
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FIGURE  II.  COORDINATE  SYSTEMS  AND  BODY  GEOMETRY 


FIGURE  12.  BEHAVIOUR  OF  THE  RECOVERY  FACTOR 
AT  THE  WALL  WITH  FV_ 
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FWJRC  13.  BEHAVIOUR  OF  THE  LOCAL  RECOVERY  FACTOR 
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FMUftC  14.  STAGNATION  POINT  HEAT  TRANSFER 
COMPARISON 
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FIGURE  15.  STAGNATION  POINT  HEAT  TRANSFER  CORRELATION 


6k 


65 


4 


SI 

©  H 
(0  O 

«S° 

*8 

M  H 


^  fl  O, 

H  O  h  I 


0) 

(h 

o 

,1 

I  ® 

:  £ 

'  H  - 


CO  O 

6  ° 


^7 


& 

® 


« 


I 

NO 

3 


fe  * 
S  ©  1 


®  © 

C  ‘ 

*H  _ 

B  ►* 
©  « 

H  «H 


2*  P  r> 


j  E  ^  .  _ 

s  S  &  S  -c 


X  3 

SI 

©  H 
©mo 

■P  ** 
CO  S  C 
P  o 
o  *H  *H 
•h  «  p 

£  (h  e 

O 


©  © 

a  a 

•  CO  o 

ego 

£  uS 
°^& 


&1 1?  §$  s  ga^sis^ic 


•  • 
a  s 
i  S  8t 

©  ®  a  o 
ho  p 
p  P.  «  ® 

£  -p  >»? 

•  o  ® 

S 

•H  >  8  | 

P  8 

95^^ 

srs 

(h  1h  O 
®  ©  >  © 
G  P  ®  H 
©  G  -H 

3  ©  Cm 

a  o  o  o 

G  U  tc,  U 

me  ex 


a  § 
u  0  *tl 

°?8 

■PH® 

8  Cm  H 

« 

Cm  H  m 

•8* 
©  H  m 
P  >  O 
-PC® 
•H  P 
Cm  © 
O  ©  B 
P  ® 
a  p  o 

•H  « 

a  c  ex 
S*H 
rH  ^ 


©  (0 


I 


§  fc,“ 

m  ex  o 

sg-° 

Cm  -p  © 
^  C  £ 
**  ©  P 


/  N 

v  I 


t 


^7-' 


A  , 

b£ 

S-.  c 
•H  ^ 

P 

O'® 

.S’-1 

m  m 
>  *H 

s  i 


£ 


B3B" 

Eh  <  fc-  >0 


l 

OS 


©  -H  I  _  m  . 

•S£5SSg|  . 

aa  c  xU  a  © 

S^rasii 

cfewH  H.r> 


CO 

H  (kt 

a 

o 

1  .. 

:  o 


t5  -r, 

5  9 
ZS 

C 


o  <  H  02  © 

E pBS  *•  - 

©  5R  f 

g|lk§£3 


i  M  re  m 

:$is* 


+»  .d 

a  $ 

n 

8  fr* 

«  -P 
J3  •£  © 

5  Id 

•ss* 

•c  h  § 
©  «  H 
p  c  p 

C  *rt 

m  3  „ 

©HO 

u  o 

(X  ©  _ 
to  2  q 


iu , 

©  a  *h  , 
o»2w/ 
ex  to 
a  ©  v 

c  >  g. 

■ .  § 

©  o 

a  s , 

p  «  x 
©OP 
o 

©  to  © 

P«-H  > 

fi  «  H 

©  s  o 

PH® 

« 

©CO 

cap 


•  ~  -  !  § 


Cm  © 

a  ' 


an© 
■H  ©  *H 

nee 

o 


s  s 

a  m 

53l 


H  -S 
«  P 

9 

■5  i  B 

P  P  H 
•HO® 
o  ®  0 
y  h  So 

a  <c  ® 

©  ®  p 

a.a>3 


9 

a  s  „ . 

S  jt* 

9  O  O  O  > 


o 

d  5 


PH  8 


P  t*>" 
SO  C  t 


•H  > 

Is. 


53  *5 

u  * 


©  H  ® 

P  ^  £} 
PC® 
•H  P 


JS8 

a  p  O 


I 


•t  ©  ® 

1  ® 
.  .  •  8^ 
tr&a 

IS 

a  P  K  h  o 
3  p  p  „ 

g**H  P 

•HP 

2  SH  . 
op  ©  j» 

G  P  ®  H  ®  C  © 
a  a  «h  ©  h 
9  9  Cm  O  H 
a  o  o  p  9  t  >* 
o  o  •  h  *t  o  H 

•g®ft§fei^ 

8  so ^ p  s § 

S  5  *8  Sgi 

a  <r  0  *H  m 

3  *  H  %4  P  © 

“  *  OP 

9  P 


/  ' 

V  / 


t 


